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5h ' Abstract. We present a large set of theoretical isochrones, whose distinctive features mostly reside on the greatly improved 

C/3 I treatment of the thermally pulsing asymptotic giant branch (TP-AGB) phase. Essentially, we have coupled the TP-AGB tracks 

Cu . described in Paper I, at their stages of pre-flash quiescent H-shell burning, with the evolutionary tracks for the previous evo- 
lutionary phases from Girardi et al. (2000). Theoretical isochrones for any intermediate value of age and metallicity are then 

#vj , derived by interpolation in the grids. We take care that the isochrones keep, to a good level of detail, the several peculiarities 

^ ■ present in these TP-AGB tracks - e.g. the cool tails of C-type stars owing to the use of proper molecular opacities as convective 

fSJ ' dredge up occurs along the TP-AGB, the bell-shaped sequences in the Hertzsprung-Russell (HR) diagram for stars with hot- 

^Nj , bottom burning, the changes of pulsation mode between fundamental and first overtone, the sudden changes of mean mass-loss 

0^ ' rates as the surface chemistry changes from M- to C-type, etc. Theoretical isochrones are then converted to about 20 different 

\l I photometric systems - including traditional ground-based systems, and those of recent major wide-field surveys such as SDSS, 

,__j ■ OGLE, DENIS, 2MASS, UKIDSS, etc., - by means of synthetic photometry applied to an updated library of stellar spectra, 

,-H ' suitably extended to include C-type stars. Finally, we correct the predicted photometry by the effect of circumstellar dust during 

f — , the mass-losing stages of the AGB evolution, which allows us to improve the results for the optical-to-infrared systems, and 

("^ ■ to simulate mid- and far-IR systems such as those of Spitzer and AKARI. We illustrate the most striking properties of these 

►-^ ' isochrones by means of basic comparisons with observational data for the Milky Way disk and the Magellanic Clouds. Access 

• 1— ( , to the data is provided both via a web repository of static tables (http://stev.oapd.inaf.it/dustyAGB87 and CDS), and 

J^ ' via an interactive web interface (http://stev.oapd.inaf.it/cmd) that provides tables for any intermediate value of age 

Vh ' and metallicity, for several photometric systems, and for different choices of dust properties. 

1. Introduction thermally-pulsing asymptotic giant branch phase (TP-AGB). 

The coverage of the TP-AGB evolution, however, has always 

Theoretical isochrones find a wide range of applications in ^^^^ f^^ ^^^^ ^^^^^^^ ^^^ physically accurate, as we will illus- 

modern astrophysics, going from the simple age-dating of star ^^^^^ ^^^^^ j^^ p^^^ modelling of this phase greatly contrasts 

clusters, to more complex problems like the derivation of star ^-^^^ -^^ importance to the integrated Ught of stellar populations 

formation histories of resolved galaxies, and, since Chariot & ^-^^^ ^^^^ ^^^^ ^ f^^ jqS y^ ^^p ^^ ^^e Hubble time. According 

Bruzual (1991), the modelling of their spectrophotometnc evo- ^^ p^^g^j ^^ ^j (^990), the TP-AGB contribution to the total 

lution across the Hubble time. luminosity of single-burst stellar populations reaches a maxi- 

In fact, many diff^erent groups are publishing large grids of ^^^^ of about 40 percent at ages from 1 to 3 Gyr, falling to 

stellar isochrones, covering wide enough a range of stellar pa- less than 10 percent at 10 Gyr. Obviously, TP-AGB stars have 

rameters (age, mass, metal content, alpha-enhancement, etc.) a great impact on the rest-frame near-infrared spectra of gala- 

to become useful in the stellar population synthesis of galaxies, xies. Moreover, TP-AGB stars account for most of the bright- 

These isochrones normally comprehend a detailed description infrai-ed (IR) objects in resolved galaxies, as clearly demon- 

of all evolutionary phases from the zero-age main sequence strated by DENIS, 2MASS, SAGE and S^MC data for the 

(ZAMS) up to either carbon ignition or to the onset of the Magellanic Clouds (Cioni et al. 1999; Nikolaev & Weinberg 

2000; Blum et al. 2006; Bolatto et al. 2007). 
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Table 1. Available stellar isochrones including the TP-AGB phase. 



Reference 



yDV and HBB*"* 



J^ 



Dust*'^* public www address for retrieval 



Bertellietal. (1994) 


NO 


solar-scaled comp. 


NO 


Bressanetal. (1998) 


NO 


solar scaled comp. 


YES 


Girardi et al. (2000) 


NO 


solar-scaled comp. 


NO 


Marigo & Girardi (2001) 


YES 


solar-scaled comp. 


NO 


Mouhcine (2002) 


YES 


solar scaled comp. 


YES 


Piovan et al. (2003) 


NO 


solar scaled comp. 


YES 


Marigo et al. (2003a) 


YES 


0-/C-rich comp. 


NO 


Cioni et al. (2006) 


YES 


0-/C-rich comp. 


NO 


Cordier et al. (2007) 


NO 


solar scaled comp. 


NO 


This work 


YES 


0-/C-rich comp. 


YES 



http://pleia di.pd. as tro.it| 
upon request 



http://pleiadi.pd.astro.it 

http://pleiadi.pd.astro.it 

upon request 

upon request 

upon request 

http://pleiadi.pd.astro.it 

http://www.oa-teramo. inaf.it/BASTI, 

http ://stev. oapd.inaf . it/cmd 



*°' Explicit inclusion or not of AGB nucleosynthesis, i.e. the third dredge-up and hot-bottom burning, in the underlying AGB models 

** Choice of low-temperature opacities affecting the predicted effective temperature, i.e. for solar-scaled or variable (O- or C-rich) chemical 

mixtures 

*'^' Inclusion or not of circumstellar dust in the computation of the emitted spectrum 



The situation regarding the TP-AGB phase in published 
sets of isochrones is the following: Whereas early isochrones 
sets do not even include the TP-AGB phase, more recent ones 
do it, but still in a very approximate way, ignoring crucial as- 
pects of the TP-AGB evolution. A summary table of avail- 
able isochrones is provided in Marigo (2007) and updated in 
Table [1] herewith. In short, BertelH et al. (1994) and Girardi et 
al. (2000; hereafter GiOO) adopt a very simple TP-AGB syn- 
thetic model that does not include the third dredge-up and 
the over-luminosity above the core mass-luminosity relation 
caused by hot bottom burning (HBB). A luminosity correction 
for stars undergoing HBB is instead included by Cordier et al. 
(2007), but in this case the third dredge up - together with its 
most remarkable consequence, i.e. the formation of C stars - 
and the HBB nucleosynthesis are completely missing. Marigo 
& Girardi (2001) present the first isochrones based on tracks 
(Marigo et al. 1999) with the third dredge up included, and 
moreover its efficiency is calibrated to reproduce the C-star lu- 
minosity functions in both Magellanic Clouds (MCU Another 
distinctive aspect of Marigo et al. (1999) tracks, is that the HBB 
efficiency, together with the model effective temperatures, are 
derived from numerical integrations of the envelope structure. 

Besides the above mentioned works, the literature abounds 
in descriptions of evolutionary population synthesis models 
which have included the complete TP-AGB evolution (e.g. 
Chariot & Bruzual 2003; Raimondo et al. 2005; Vazquez & 
Leitherer 2005; only to mention a few). However, in all these 
cases, critical processes of TP-AGB evolution - again the third 
dredge-up, HBB, or superwind-like mass loss - are missing 
in the models. In some cases (e.g. Vazquez & Leitherer 2005; 
Chariot & Bruzual 2003), the tracks for the TP-AGB and previ- 
ous evolutionary phases come from different sources; this im- 
pUes that basic stellar properties like the core mass and enve- 
lope chemical composition are discontinuous along the com- 



posite tracks. As demonstrated by Marigo & Girardi (2001), 
discontinuities in these quantities also produce an energetic in- 
consistency, in the sense that they violate the strict coupling 
between stellar integrated light and chemical yields implied by 
the condition of energy conservation. 

A few authors have computed isochrones aimed at the mod- 
elling of spectral evolution of galaxies which consider the effect 
of circumstellar dust around AGB stars (Bressan et al. 1998; 
Mouhcine 2002; Piovan et al. 2003), but without improving the 
treatment of the underlying TP-AGB evolution, with respect to 
those already mentioned. 

In all the previously-mentioned cases, moreover, the effect 
of variable opacities along the TP-AGB evolution is completely 
ignored: The basic stellar properties, and in particular the ef- 
fective temperature Jeff, are derived using opacities valid for 
scaled-solar compositions, which are highly inadequate to de- 
scribe the chemically evolved envelopes of TP-AGB stars. The 
proper consideration of variable opacities is now known to be 
crucial (see Marigo 2002; Marigo et al. 2003a), especially for 
stars undergoing the transition to the C-rich phase - i.e. com- 
prising, for MC metallicities, all stars with initial masses be- 
tween about 1.5 and 3.5 Mq, and all isochrones with ages be- 
tween about 0.5 and 5 Gyr. This situation is partially reme- 
died by Cioni et al. (2006), who has made pubUcly available 
the first isochrone sets in which variable molecular opacities 
were considered along the TP-AGB evolutiorB In fact, the C- 
star phase in these isochrones is characterized by much lower 
temperatures (reaching T^ff as low as 2300 K) than former 
ones. Those specific TP-AGB models, however, are previous 
to the extensive calibration of dredge up efficiency performed 
in Marigo & Girardi (2007, hereafter Paper I). In this latter 
work, the dredge up parameters are tuned so that the models 
reproduce the observed the C-star luminosity functions in the 



' Many of the improvements in Marigo et al. ( 1 999) TP-AGB tracks 
were already introduced by Groenewegen & de Jong (1993, 1994abc), 
who computed extended sets of tracks and performed population syn- 
thesis of resolved AGB stellar populations. Groenewegen & de Jong, 
however, do not provide theoretical isochrones from their tracks. 



- Marigo et al. (2003a) also computed isochrones containing essen- 
tially the same kind of TP-AGB tracks as in Cioni et al. (2006), but for 
a more limited range of metallicities, and for exploratory choices of 
model parameters. These isochrones were provided a few times upon 
request, and were superseded first by Cioni et al.'s and now by our 
new isochrones. 
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Magellanic Clouds, and the TP-AGB lifetimes of M and C stars 
in Magellanic Cloud star clusters (Girardi & Marigo 2007a). 

The goal of the present paper is to remedy this highly un- 
satisfactory situation, providing the first set of isochrones in 
which the TP-AGB evolution is treated with a satisfactory level 
of completeness and detail - i.e. considering the crucial effects 
of third dredge up, hot bottom burning, and variable molecular 
opacities - and calibrated at the same time - i.e. able to repro- 
duce the basic observables ofAGB stars in the MCs. Moreover, 
we wish to provide these results in many photometric systems 
and including the reprocessing of radiation by circumstellar 
dust in mass-losing stars, which is crucial for the modelling of 
the infrared photometry. 

The properties of our stellar evolutionary tracks are briefly 
recalled in Sect. |2] which discusses the basic aspects of con- 
structing isochrones using the new TP-AGB tracks of Paper I 
- including new variables like the C/O ratio, pulsation periods, 
pulsation modes, etc. Section[3]deals with the transformation of 
the isochrones to many photometric systems, which require up- 
dating and extending the libraries of stellar spectra previously 
in use. The inclusion of dust reprocessing in the predicted pho- 
tometry is illustrated in Sect.|4] A preliminary comparison with 
observational data is done in Sect. |5] Section |6] describes how 
to access the data tables - containing isochrones, bolometric 
corrections, dust attenuation curves, integrated magnitudes of 
single-burst populations, etc. - and summarizes both the main 
novelties and caveats of our new isochrones. 

2. From evolutionary tracks to isochrones 

2.1. Putting tracks together 

The evolutionary tracks used in this work have already been 
thoroughly presented in previous papers. The evolution from 
the ZAMS up to the first thermal pulse on the AGE is taken 
from GiOO for the initial metaflicities Z = 0.0004, 0.001, 0.004, 
0.008, 0.019 (the solar value), 0.03, and from Girardi (2001, 
unpubHshed; see http://pleiadi.oapd.inaf.it) for Z = 
0.0001. The range of initial masses of these tracks goes from 
Mi - 0.15 to 7 Mq, out of which afl models between Mi ~ 0.5 
and 5 Mq are considered to undergo the TP-AGB phase. 

The TP-AGB phase is then computed via a synthetic code 
as described in Paper I, starting from the physical conditions at 
the first thermal pulse, and up to the complete loss of the en- 
velope via stellar winds. The model takes into account many 
critical aspects of the TP-AGB evolution, like the third dredge- 
up events, HBB over-luminosity and nucleosynthesis, pulse 
cycle luminosity and Teif variations, opacity variations driven 
by the changing chemical composition in the envelope, long- 
period variability (LPV) in both first overtone and fundamen- 
tal modes, and their modulation of the mass-loss rates M via 
changes in the pulsation period P. 

Especially relevant for the present paper are the formalisms 
adopted for the mass loss: besides affecting the AGB lifetimes, 
high mass-loss rates determine the formation of circumstel- 
lar dust and hence the obscuration of AGB stars in the opti- 
cal and their brightening in the mid- and far-IR (see Sect. Hji. 
As detailed in Paper I, for O-rich AGB stars we use the re- 



sults of dynamical atmospheres including dust for long-period 
variables calculated by Bowen (1988) and Bowen & Willson 
(1991), taking into account whether pulsation occurs on the 
first overtone or on the fundamental mode. For C-rich stars, 
we assume two diff'erent regimes of mass loss: initially stars 
follow the semi-empirical relation from Schroder et al. (2003), 
then switching, above a critical luminosity, to the M values de- 
rived from the pulsating dust-driven wind models by Winters 
et al. (2000, 2003). In both cases, mass-loss rates depend on 
the pulsation period, which is predicted along the evolutionary 
tracks according to the period-mass-radius relations derived by 
Ostlie & Cox (1986), Fox & Wood (1982), Wood et al. (1983). 

The most massive TP-AGB tracks (M; > 3.0 Mo) of Paper I 
have been replaced by new calculations, accounting for a bet- 
ter description of the transition luminosity between the first- 
overtone and the fundamental-mode pulsation. This aspect is 
thoroughly described in the Appendix. 

The synthetic TP-AGB evolution of Paper I has been fol- 
lowed up to the complete ejection of the stellar envelopes, 
hence including a good fraction of the post- AGB phase prior 
to the planetary nebula regime, i.e. the fast terminal stages dur- 
ing which the star evolves at almost constant luminosity, while 
reft rapidly increases due to the combined effect of the out- 
ward displacement of the H-burning shell, and the stripping of 
the residual envelope by stellar winds. These final parts of the 
tracks may not be reliable since some relations used in the syn- 
thetic TP-AGB code become inaccurate. Moreover, the mass- 
loss rates are critical in determining the duration of this post- 
AGB transition phase, but they are still poorly known. 

Anyway, this transition phase turns out to be very short- 
living ($ 10"* yr for low-mass stars, and even shorter for 
intermediate-mass ones) if compared to the previous TP-AGB 
evolution. Therefore, although this post-AGB phase is mostly 
included in the new isochrones, its overall impact on them is 
almost negligible. 

Another important detail is that before building the 
isochrones, TP-AGB tracks are "cleaned" of their luminos- 
ity and Teff variations caused by thermal pulse cycles. Just 
the stages of quiescent H-shell burning previous to He-shell 
flashes, along with the stage of coolest T^ff, are kept for the 
isochrone interpolation. This is a necessary approximation 
since the large luminosity and T^ff variations between confin- 
ing tracks would lead to wild variations in the interpolated 
isochrones. Anyway, the detailed evolution of luminosity and 
Teft over the thermal pulse cycles can be easily introduced in 
the tabulated isochrones a posteriori, since the shape of ther- 
mal pulse cycles is a known function of other quantities that are 
also interpolated and stored along the isochrones, like the core 
mass and surface chemical composition (see e.g. Wagenhuber 
& Groenewegen 1998). 

The GiOO and Paper I sets of evolutionary tracks diff'er in 
their source of low-temperature (T ^ 10"* K) opacities: while 
GiOO uses tables for scaled-solar chemical compositions taken 
from Alexander & Ferguson (1994), in Paper I the opacities 
are computed "on-the-fly" using the proper chemical mixture 
of the evolving TP-AGB envelopes, by means of the meth- 
ods and approximations described in Marigo (2002). There 
is no dramatic change in opacities, however, in passing from 
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the scaled-solar tables of Alexander & Ferguson (1994) to 
Marigo's (2002) formalism, provided that the passage occurs at 
C/O ratios much smaller than 1, as indeed is the case at the first 
thermal pulse of all TP-AGB tracks. The molecular opacities 
depart from those in Alexander & Ferguson (1994) only later 
in the TP-AGB evolution, when the third dredge up sets in and 
the C/O ratio increases above 1. Therefore, the concatenation 
of both sets of tracks at the stage of first thermal pulse produces 
just modest discontinuities in the tracks, i.e. small jumps in Teff 
that are caused by the small differences in opacities for O-rich 
mixture^ For other quantities, like the luminosity, core mass, 
and surface chemical composition, the inspection of isochrones 
does not reveal any discontinuity. 

Although the reader may feel uncomfortable with the small 
discontinuities in Tgff at the stage of first thermal pulse in the 
isochrones, we remark that they are anyway dramatically less 
important than any discontinuity that one may find when amal- 
gamating tracks from completely different sources. The simple 
fact that the core mass and envelope chemical composition vary 
in a smooth and continuous way, for instance, ensures that the 
total emitted light derived from a given set of isochrones will 
be fully consistent with the chemical yields and total remnant 
masses derived from the same sets of tracks. As illustrated by 
Marigo & Girardi (2001), this is an important, although largely 
ignored, consistency requirement to evolutionary population 
synthesis models of galaxies. 

2.2. Resulting isochrones 

Theoretical isochrones are built via simple linear interpola- 
tions inside the grids of evolutionary tracks, performed be- 
tween "equivalent evolutionary points". The final isochrone, 
for a given age t, is made of a sample of such points conve- 
niently distributed on the HR diagram. For each point, several 
stellar parameters are stored. In fact, the newest TP-AGB tracks 
take into account various aspects of evolution that have been 
ignored while making the GiOO isochrones, but that are essen- 
tial components of the TP-AGB models of Paper I, namely: 
surface chemical composition (especially the C/O ratio), stel- 
lar core and envelope mass, mass-loss rates, pulsation modes 
and periods. These parameters can, for instance, allow an easy 
identification of C-type and dust-obscured stars, and are essen- 
tial to the simulation of variations in stellar properties due to 
thermal pulse cycles and long-period variability. Therefore, our 
isochrone-making routines have been extended to take explicit 
care of these properties. 

When an isochrone is built out of our tabulated tracks, 
mass loss and LPV pulsation are taken into detailed account 
only after the transition from the early- AGB to the TP-AGB 
regime. Variables such as the mass-loss rates, pulsation peri- 
ods, and circumstellar dust emission, pass from "non-defined" 
values to well defined ones, corresponding to the occurrence 
of the first thermal pulse. In order to avoid discontinuities in 
these variables, we have a posteriori evaluated the mass-loss 



^ These jumps in Tcti can be appreciated in Fig. \T\ at the points 
where the TP-AGB phase begins (i.e. in passing from green to blue 
lines). 



rates and pulsation properties of stars in earlier phases, even 
if these have been computed at constant mass. The procedure 
is straightforward: for each star in the pre-TP-AGB section of 
the isochrones, the basic stellar properties (L, Tetf, Mi) are used 
to determine pulsation periods and mass-loss rates in the same 
way as for O-rich AGB stars. Notice that in most cases mass 
loss prior to the TP-AGB phase is expected to be very modest 
and not to have evolutionary implications, so that this proce- 
dure is well justified. The only exception is for low-mass mod- 
els close to the tip of the red giant branch (TRGB), where the 
use of Reimers (1975) mass-loss formula with a multiplicative 
factor of 77 = 0.4, leads to the loss of AM^'^^ ^ 0.2 Mq for 
low-mass starO This amount of mass lost is necessary to ex- 
plain the HB morphologies observed in galactic globular clus- 
ters (see Renzini & Fusi Pecci 1988). The use of the Reimers 
mass-loss formula for the RGB shall be seen as just a "con- 
servative choice" since, although alternative formulas do ex- 
ist (see Catelan 2000, 2007), at the present time none can be 
clearly favoured over the others. We take the RGB mass loss 
into consideration during the construction of stellar isochrones 
by simply jumping from a TRGB track of given mass M to a 
ZAHB track of mass M - AM^'^^. The procedure is somewhat 
detailed in BertelH et al. (1994) and in GiOO. 

Extended sets of isochrones are depicted in the HR dia- 
grams of Fig. [T] Since the main novelty of the present paper 
is in the detailed description of the TP-AGB phase, the HR di- 
agrams are zoomed on this phase. In other parts of the HR, 
isochrones are simply identical to those previously computed 
by GiOO. 

It is evident from the figure that our new isochrones reflect 
all of the peculiar behaviours of the TP-AGB tracks they derive 
from (see figure 20 in Paper I). In particular, one can notice the 
following: 

At moderate-to-high metallicities, i.e. for Z > 0.004 : In the 
youngest isochrones, the combination of HBB and high mass- 
loss rates gives origin to bell-shaped TP-AGB curves. The lu- 
minosity first rapidly grows above the core-mass luminosity re- 
lation because of the occurrence of HBB, then higher mass-loss 
rates are attained, drastically reducing the stellar envelope mass 
so that HBB and the related over-luminosity are eventually ex- 
tinguished. 

These stars usually have C/0< 1, although a late transi- 
tion to the C-rich phase may take place due to final dredge up 
events when HBB has already ceased (see Frost et al. 1998) in 
agreement with observations (e.g. van Loon et al. 1998, 1999). 

The C-star phase is completely developed in isochrones of 
intermediate ages, and it happens at very low Teff (reaching 
reft - 2500 K), in a sort of "red tail" which is neatly separated 
from the sequences of O-rich stars. This separation naturaUy 
results from the very different low-temperature opacities of C- 
rich envelopes (Marigo 2002). 

In the oldest isochrones, dredge up is not operating and the 
TP-AGB sections of the isochrones simply align above the pre- 
vious early-AGB phase. They steadily increase in luminosity 



■* The behaviour of AM''°^ as a function of age can be appreciated 
in figure 6 of Girardi (1999), top panel. 
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Fig. 1. Isochrones in the H-R diagram , zooming on the region of the AGB, for 
the 7 metaUicity values of the original tracks. The isochrones are plotted with ages 
log(f/yr) - 7.8 to 10.2 (the youngest are at the top/left) with a spacing of 0.1 dex 
in log t. In the electronic version of the paper, the pre-TP-AGB phases are drawn in 
green, whereas on the TP-AGB the surface O-rich and C-rich configurations are dis- 
tinguished with blue and red lines, respectively. Along the TP-AGB, the luminosity 
and effective temperature correspond to the quiescent values just before the predicted 
occurrence of thermal pulses. 



J 7 3.6 3.5 3.4 

log T^„ 



and slightly bend to the left at their upper end, as a result of the 
progressive thinning of stellar envelopes. 



At low metallicities, i.e. for Z < 0.004: Some of the previ- 
ous trends change. One important difference is that, despite 
the high efficiency of HBB, a C-star phase is present even at 
younger ages. This is caused by the activation of the ON cycle 
which converts O into N, so that the C/O ratio is predicted to in- 
crease above unity (Ventura et al. 2002; and Paper I). Compared 
to the tracks of higher metallicities, the C-star sequences start 
at higher effective temperatures, as expected from the lower 
molecular concentrations in their atmospheres. The oldest O- 
rich sequences evolve to larger luminosities while bending to 



the left at their high ends, just as in the tracks for higher metal- 
licities. This feature is essentially related to the outward dis- 
placement of the H-burning shell, which is found to be more 
efficient than the stellar winds in reducing the envelope mass 
of low-mass AGB stars {M\ S, 1. Mq). 

A more detailed discussion of the above-mentioned be- 
haviours is presented in Paper I. 

In order to better illustrate the peculiarities of the TP- 
AGB models presently available. Fig. |2] presents a series of 
Z = 0.008 isochrones derived from GiOO tracks connected 
with different extensions of the TP-AGB phase. In the first one 
(GiOO) the TP-AGB phase is very simplified and just extends 
the previous evolution straight above the RGB and early- AGB 
phases, without any C-star sequence and without any signa- 



Marigo et al.: Optical to far-IR isochrones 




3.6 3.5 

log T,„ 




Fig. 2. The same as in Fig. [T] but for a sequence that illustrates the evolution of the TP-AGB phase since the release of GiOO 
isochrones. From left to right, isochrones are from GiOO, Marigo & Girardi (2001), and Cioni et al. (2006). This sequence would 
be completed with the third panel of Fig. [T] 



ture of HBB. Then, the Marigo & Girardi (2001) isochrones 
do introduce dredge-up and HBB in the TP-AGB, and present 
a well-defined C-star region, which is still aligned on top of 
the RGB and early-AGB, as expected from models using opac- 
ities for scaled-solar compositions. Cioni et al. (2006) intro- 
duce variable molecular opacities and then the C-star sequence 
clearly changes its shape on the HR diagram (as explained by 
Marigo 2002). Finally, the sequence of Fig. |2] could be com- 
pleted with the third panel of Fig. [1] which represents the 
present state of the Z = 0.008 isochrones, differing from the 
Cioni et al. (2006) in that the C-star luminosities and lifetimes 
are now calibrated (see Paper I). 

3. Transformations to several photometric 
systems 

The new isochrones are transformed to several photometric 
systems, following the formalism described in Girardi et al. 
(2002) to derive bolometric corrections from a library of stel- 
lar spectra. Although this argument will be more thoroughly 
discussed in a forthcoming paper (Girardi et al. 2007, in prepa- 
ration), it is worth recalling the main steps we follow to deal 
with the novelties of the present release. 

3.1. Bolometric corrections and extinction coefficients 

Given the spectral flux at the stellar surface, F^, bolometric 
corrections for any set of filter transmission curves S a are given 
by (see Girardi et al. 2002 for details) 



face. Once BCs , are computed, stellar absolute magnitudes fol- 
low from 



BCs 



Mboi,o - 2.5 log [47r(lOpc)Vboi/L0] 

•J A] 



(1) 



-1-2.5 log 



it^f'SAdA 



where /° represents a reference spectrum (at the Earth) that 
produces a known apparent magnitude m" , and Fboi = 

J. F,idA - crT^g is the total emerging flux at the stellar sur- 



Ms, = Mboi - BCs, , 
where 

Mboi = Mboi.G - 2.5 log(L/L0) 

= Mboi,0 - 2.5 log(4;r/?Vboi/L0) 



(2) 



(3) 



As in Girardi et al. (2002), we adopt Mboi,© = 4.77, and Lq - 
3.844 X 10^"* erg s"' (Bahcall et al. 1995). 

Notice that Eq. ([T]i uses photon count integration instead 
of energy integration, as adequate to most modern photometric 
systems which use photon-counting devices instead of energy- 
amplifiers. 

The equation also includes the interstellar extinction curve 
Aa, that we set to A^ = 0. For the optical and near-IR, we pro- 
vide the extinction coefficients A^/Av as computed for a yel- 
low dwarf of solar metallicity (the Sun) using the Cardelli et 
al. (1989) Rv - 3.1 extinction curve. The small variations of 
such coefficients with spectral type and extinction (see Grebel 
& Roberts 1995) will be discussed in a separate paper (Girardi 
et al., in prep.) which will detail the transformations to WFPC2 
and ACS systems, and describe a way to consistently include 
the interstellar extinction into the isochrones. 

3.2. Assembling a library of spectral fluxes 

Girardi et al. (2002) describes the assembly of a large library 
of spectral fluxes, covering wide ranges of initial metallicities 
([M/H], from -2.5 to -1-0.5), Tes (from 600 to 50000 K), and 
log^ (from -2 to 5JS This grid is wide enough to include 
most of the spectral types that constitute the bulk of observed 
samples. It is composed of ATLAS9 "NOVER" spectra from 
CastelU et al. (1997; see also BesseU et al. 1998) for most of 
the stellar types, Allard et al. (2000; see also Chabrier et al. 



In this paper, the surface gravity g is expressed in c.g.s. units. 
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2000) for M, L and T dwarffl Fluks et al. (1994) empirical 
spectra for M giants, and finally pure black-body spectra for 
stars exceeding 50 000 K. 

Since then this library been revised. We have replaced the 
ATLAS9 spectra by the newest "ODFNEW" ones from Castelli 
& Kurucz (2003), which incorporate several corrections to the 
atomic and molecular data, especially regarding the molecu- 
lar line Usts of CN, OH, SiO, H2O and TiO, and the inclusion 
of HI-H-H and HI-Hh- quasi-molecular absorption in the UV. 
These new models are considerably better than previous ones 
in the UV region of the spectrum, and over most of the wave- 
length range for Teff between ~ 4500 and 3500 K (see figs. 1 to 
3 in Castelli & Kurucz 2003). 

Moreover, the library has been extended with the Loidl et 
al. (2001) synthetic C-star spectra for 3600 k (T^g/K) > 2600, 
derived from hydrostatic model atmospheres computed with 
MARCS (see also Hofner et al. 2003 for more details). For 
C-stars cooler than 2600 K, we simply adopt the BC values 
derived from 2600 K spectra. These Loidl et al. (2001) spectra 
are computed for C/0= 1 . 1 for all Teff, and for both C/0= 1 . 1 
and 1.4 for Teff < 3200 K. We have verified that bolometric 
corrections for all broad bands redder than V are quite insensi- 
tive to the C/0 ratio, varying by less than 0. 1 mag between the 
C/0=1.1 and 1.4 cases. These differences become S 0.01 mag 
in the 7-band. Just for the bluest optical pass-bands the choice 
between these two C/O values causes significant differences the 
bolometric corrections, amounting for instance to ^ 0.4 mag in 
B. Anyway the Teff < 3200 K C-stars are much fainter (more 
than 4 mag in a Vega system) at these blue wavelengths than 
at the /-band, and therefore these cases are of little practical 
interest. Hence, we adopt by default the C/0=1.1 C-star spec- 
tra, for all Teff. A more elaborated approach will be followed 
after spectra for a wider range of C/O values become available 
(Aringer et al., in preparation). 

Finally, as part of the TRILEGAL project (Girardi et al. 
2005), we have included spectra for DA white dwarfs from 
Finley et al. (1997) and Homeier et al. (1998) for T^ff between 
100 000 and 5 000 K, and log g between 7 and 9. 

The CastelU & Kurucz (2003) and Allard et al. (2000) tabu- 
lated spectra extend well into the far-IR, until 160 and 971 /im 
respectively. In the other cases, the available spectra were tab- 
ulated just up to the mid-IR, namely until 12.5 jum for Fluks 
et al. (1994) and 2.5 fim for Loidl et al. (2001). In all cases 
we have extended the flux to the very far-IR using a scaled 
Rayleigh-Jeans law starting from the point of maximum tabu- 
lated A. This is clearly a crude approximation, van der Bliek 
et al. (1996), for instance, find that black-body fluxes in the 
far-IR present discrepancies of the order of 20 per cent with 
respect to the solar flux. Decin & Eriksson (2007) recently dis- 
cussed the status of the predicted spectral energy distributions 
in the IR, to whom we refer for a detailed discussion of this 
topic. It turns out that our predicted BC for the mid- and far-IR 
will be uncertain by up to a few tenths of a magnitude for the 



* Although brown dwarfs are not the subject of this paper, the tables 
of bolometric corrections we are going to distribute do include these 
objects. More details can be found in the original paper by Allard et 
al. (and to a smaller extent also Girardi et al. 2002). 



coolest M and C stars. However, we expect that the errors in 
colours will be somewhat smaller, and they will become irrele- 
vant in stars with significant mass-loss, which are exactly those 
with the most prominent fluxes in the mid- and far-IR. In fact, 
in these cases the emission at mid- and far-IR wavelengths is 
dominated by the circumstellar dust shells and not by the stellar 
photosphere (see Sect. |4|i. 

3.3. Filter sets so far considered 

Table |2] presents all the photometric systems we have con- 
sidered so far. They include several traditional ground-based 
photometric systems (Johnson-Cousins-Glass, Stromgren- 
Crawford, Washington, BATC), systems belonging to some 
influential wide-field surveys (DENIS, 2MASS, OGLE-II, 
SDSS) and cameras (WFI, UKIDSS), and then to some of 
the most successful satellite-based cameras of the recent past 
(Tycho, ACS) or still in operation (WFPC2, NICMOS, Spitzer, 
AKARI). Other systems are presently being added to the list. 

Note that the transformations to some of these photometric 
systems are described in other papers (e.g. Girardi et al. 2007 
for WFPC2 and ACS; Vanhollebeke et al. 2007 for OGLE-II); 
they are listed in Table|2]just for the sake of completeness and 
no further detail on these systems will be provided here. All 
systems in the Table are already available at the web interface 
(see Sect. 16. 11 1. 

4. Reprocessing by circumstellar dust 

Circumstellar dust reprocesses a substantial fraction of the pho- 
tons emitted by the coolest evolved stars. Since, on one hand, 
the efficiency of dust extinction decreases quickly above ~ Ifim 
and, on the other hand, dust does not survive at temperatures 
above 1000-2000 K, the general effect of the presence of dust 
in an astrophysical source is to shift part of the intrinsic emitted 
power from the optical-ultraviolet to the mid- and far-infrared 
region. 

A serious difficulty in taking into account in detail the ef- 
fects of dust is that the optical properties (i.e. the efficiency 
of absorption and scattering as a function of wavelength) dif- 
fer significantly in different environments. This is particularly 
true for the envelopes of AGB stars, which are believed to be 
one of the two main producers of interstellar grains, the other 
being the supernovae blast waves (e.g. Dwek 1998; Todini & 
Ferrara 2001; Morgan & Edmunds 2003; Calura et al. 2007). 
As for the dust composition, the general wisdom is that the en- 
velopes of C-type stars are composed by carbonaceous grains 
while those of M-type stars are dominated by silicates, however 
several recent detailed studies, mainly based on Spitzer data 
(e.g. Lebzelter et al. 2006; Verhoelst et al. 2006) and in more 
detailed modelling (e.g. Ireland & Scholz 2006; Woitke 2006; 
Ferrarotti & Gail 2006) demonstrated that the situation may be 
more complex. Moreover, the grain size distribution, which to- 
gether with the composition determines the optical properties, 
is expected to be significantly different from that infeiTed in the 
general interstellar medium. 

A further big difficulty in the understanding of AGB en- 
velopes is the occurrence of strong mass loss and the possi- 



Table 2. Photometric systems for which our new isochrones are made available (updated to November 2007). 



Photometric system Filter names 



O 



o 

•c 



Central 
wavelength range 



kind of 
zero-points 



references 
for filters 



additional information 



Johnson 


UBV 


3642 - 5502 A 


Vega*"' 


Maiz-Apellaniz (2006) 


Bessell & Brett 


UxBxBV(Rr)c 
JHKLL'M 


3627 A- 4.72 yum 


Vega 


Bessell (1990) 
Bessell & Brett (1988) 


UBVRIJHK'^^ 


UBV(RI)cJHK 


3642 A- 2.15 yum 


Vega 


Maiz-Apellaniz (2006) 

Bessell & Brett (1988), Bessell (1990) 


Stromgren-Crawford 


uvby 


3474 - 5479 A 


Vega'"' 


Maiz-Apellaniz (2006) 
Moro & Munari (2000) 


Tycho 


VjBj 


4250 - 5303 A 


Vega'"' 


Maiz-Apellaniz (2006) 


SDSS 


ugriz 


3587 - 8934 A 


ABf-^' 


Fukugitaetal. (1996) 


OGLE-II 


UBVI 


3695 - 8321 A 


Vega 


OGLE website 


ESO/WFI 


wide filters 


3469 - 8590 A 


Vega 


ESO/WFI website 


Washington 


CMT1T2 


3982 - 8078 A 


Vega 


Geisler(1996) 


DENIS 


UK, 


7927 A- 2.15 yum 


Vega 


Fouque et al. (2000) 


2MASS 


JHK, 


1.23-2.16//m 


Vega(=> 


Cohen et al. (2003) 


UKIDSS 


ZYJHK, 


8823 A- 2.20 yum 


Vega 


Hewett et al. (2006) 


BATC 


all 15 filters 


3383 - 9722 A 


AB 


Yan et al. (2000) 


HST/WFPC2 


all wide filters 


1700 - 9092 A 


Vega+ST+AB 


Holtzmanetal. (1995) 


HST/ACS 


all WFC+HRC 


2645 -9130 A 


Vega 


Sirianni et al. (2005) 


HST/NICMOS 


7//A'-equivalent 


1.11- 2.04 //m 


Vega+ST+AB 


D. Figer, priv. comm. 


Spitzer 


all MIPS+IRAC 


3.52- 154 //m 


Vega,AB<* 


Spitzer website 


AKARI 


all filters 


2.32- 160 //m 


Vega,AB<* 


AKARI website 



supersedes Girardi et al. (2002) 



(f) 



details in Girardi et al. (2004) 
details in Vanhollebeke et al. (2007) 
supersedes Girardi et al. (2002) 
details in Girardi et al. (2002) 
details in Cioni et al. (2006) 
supersedes Bonatto et al. (2004) 
supersedes files distributed since 2006 
supersedes files distributed since 2006 
details in Girardi et al. (2007) 
details in Girardi et al. (2007) 
details in Girardi et al. (2002) 
details in Groenewegen (2006) 
details in Groenewegen (2006) 



*"' Zero-points are also taken from Maiz-Apellaniz (2006). 

*' Combines UBV from Maiz-Apellaniz (2006) with (RI)cJHK from Bessell & Brett (1988) and Bessell (1990). This combination is provided as being the most useful for the analysis of 

ground-based data. Supersedes Girardi et al. (2002). 

*'^* We adopt the traditional definition of magnitudes, using Pogson's formula, and not Lupton et al.'s (1999) one; see Girardi et al. (2004) for the motivations of this choice. 

*''' A Vegamag zero-point reproduces the photometry quite accurately for most IR systems (see Reach et al. 2005). AB magnitudes, however, may be convenient because then the magnitudes can 

be easily converted to the flux in milli-Jansky by means of a single formula, i.e. /, (mjy) = -0.4 toab (mag) - 6.56. 

*°* We use zero-points determined by Maiz Apellaniz (2007) using Bohlin's (2007) Vega spectrum. As a consequence, the 2MASS zero-point offsets measured by Cohen et al. (2003) do not 

apply 

® For the II/3 filters, we simply assume/? = H/Jwide - H/Jnarrow = 2.873 for Vega (cf. Hauck & Mermilliod 1998). The zero-point is arbitrarily set as II/3„i(ie(Vega) = 0. 
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ble relation with the type of dust that can condensate in their 
outflows. The massive winds in AGB stars are commonly as- 
cribed to the coupled eff'ect between large amplitude pulsations 
and radiation pressure on dust grains. According to detailed dy- 
namical studies, in C-rich stars the radiation pressure on dust 
can in fact be strong enough to explain the observed outflows 
(see e.g. Winters et al. 2000; Hofner et al. 2003). Instead in 
O-rich stars the possible type of grains are either not stable or 
opaque enough near ~ l//m to explain the outflows (e.g. Ireland 
& Scholz 2006; Woitke 2006). 

Furthermore, recent spectroscopic analyses of oxygen-rich 
AGB stars based on ISO data for the LMC (Dijkstra et al. 2005) 
and the Galactic Bulge (Blommaert et al. 2006) have shown a 
clear correlation between the dust composition and the mass- 
loss rates. For instance, the dust features displayed by oxygen- 
rich stars are mainly ascribed to AI2O3 for low M, while go- 
ing to higher M they are interpreted with an increasing amount 
of silicates. These findings are consistent with the latest re- 
sults from Spitzer data of AGB stars belonging to the globu- 
lar cluster 47 Tuc (Lebzelter et al. 2006). The dust mineralogy 
actually changes as the star evolves at increasing luminosities 
and pulsation periods. On the other hand, for C stars unveiling 
the dust-condensation sequence is more problematic due to its 
complex dependence on a number of factors (Speck et al. 2006; 
Thompson et al. 2006). 

What emerges is that, despite the conspicuous progress in 
both observational and theoretical grounds, the issue of dust 
formation and evolution in AGB stars is still aff'ected by large 
uncertainties. For these reasons in this paper we will explore 
just a few possible kinds of circumstellar dust, essentially fol- 
lowing the prescriptions by Groenewegen (2006, G06) and 
Bressan et al. (1998, B98). Alternative and more complex de- 
scriptions will be examined in subsequent papers. 

4.1. Method 

This section outlines the adopted procedure to include the effect 
of dust in our isochrones. 

In general, for a given star of (L, reft), losing mass at a 
rate M, with dust expansion velocity Vexp and dust-to-gas ratio 
"P, radiation transfer (RT) calculations provide the dust optical 
depth T and its corresponding attenuation/emission spectrum. 
The optical depth t scales with the other quantities via the ap- 
proximated relation 



T(ljum)^AdM*I'v;^'pL-"-^ 



(4) 



(see G06 and references therein) where Aj is a quantity that is 
set mainly by the dust composition. 

For practical convenience we have pre-computed a number 
of reference RT models which are parametrized as a function 
of a basic set of parameters, namely: T(l//m), Tf^g, L, dust con- 
densation temperature T^, and for different dust mixtures, as 
specified in Sects. 14.4. l1 and |4.4.2| for G06 and B98 models, re- 
spectively. Then, each star along the isochrone is assigned its 
proper dusty envelope by using Eq. (|4]l with the corresponding 
values of M, ^, Vexp, L, Teff, and Aj. 

Each determination of T(lyum) corresponds to specific dust 
attenuation and emission functions, which in turn fix the ratio 



between the flux derived from RT calculations and the dust- 
free intrinsic flux, i.e. F^IF^q. Finally, bolometric corrections 
are derived, i.e. ABC^ - -2.5 log(F^/F^,o)- 

In summary, for each AGB star along the isochrone, the 
general scheme proceeds through the following steps: 

- We first distinguish between O- and C-rich envelopes, us- 
ing the photospheric C/O ratio. The dust composition is se- 
lected among those for which RT calculations have been 
performed. 

- Vexp and *}* are computed as a function of L, Teft, C/O and M, 
using the formalism described in Sects. l4r2l and l4.3l below. 

- The proper t is calculated via Eq. (|4]i. 

- The corresponding ABC^ relation is interpolated at the 
effective wavelengths A^s of all filters considered in our 
isochrones. The resulting ABC^^^ values are added to the 
absolute magnitudes computed previously via Eq. ([T). 



4.2. Expansion velocity 

The wind terminal velocity or expansion velocity, Vexp, is calcu- 
lated according to the formalism presented by Elitzur & Ivezic 
(2001), to which the reader should refer for the details of the 
underlying theoretical analysis. In brief, following Elitzur & 
Ivezic (2001) the dust-wind problem is solved by means of a 
suitable scaling approach, so as it is possible to express Vgxp as 
a function of mass-loss rate, stellar luminosity and various dust 
parameters (dust-to-gas ratio, chemical composition, conden- 
sation temperature, etc.) in the form 

,4/3 N -1/2 



{AM-6} 



1/3 



1 +B- 



M" 



U 



(5) 



where v'exp is expressed in km s '; M_6 is the mass-loss rate 
in units of IO^^Mq yr"'; L4 is the stellar luminosity in units of 
104 Lo. 

The dependence of Vexp on the dust-to-gas ratio is included 
in the definition of the parameters A and B (through the quantity 
0-22, see below): 



A^3.08lO'T^,Q.crl,Xo' 



B = 



2.28 



eyvr 



^3/4 1/2^ 
(Jv ^11 ''cS 



-4/3 



(6) 



(7) 



which in turn depend on several quantities quoted here below. 
The condensation temperature in units of 10^ K is denoted with 
Tea. The parameter 2* is the Planck average of the efficiency 
coefficient for radiation pressure, evaluated at the effective tem- 
perature, while Qy is the efficiency factor for absorption at vi- 
sual. The quantity ;t'o is the ratio: 

2p(7"eff) ,„, 

where 2p(r) is the Planck average of the absorption efficiency 
at the temperature T . 

We define 0-22 - crgas/(10"^^ cm^), where (t„s& (in cm^) is 
the dust cross-sectional area per gas particle at condensation 

2 "dust fr.\ 

CTgas = na . (9) 
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In the above equation a is the mean size (in cm) of the dust 
grains, njust and Wgas are the number densities (in cm"-') of dust 
and gas particles, respectively. Similarly to Elitzur & Ivezic 
(2001), here it is assumed that all dust grains have the same 
size. 

The dust-to-gas mass ratio is defined as 



"¥ 



Pdust 
Pgas 



(10) 



where pdust and pg^s correspond to the density (in g cm"-') of 
the matter in the form of dust and gas, respectively. 

Assuming that the dust grains have a specific density pgrain 
(in g cm"-') the explicit form of the dust density reads 



Pdust - :z^a Pgiain "dust 



(11) 



With the approximation that all the gas is made up of hydrogen 
and helium, with abundances (in mass fraction) Xu and Xhs 
respectively, then 



Pgas — ^gas'^H^gas 



(12) 



where Agas - 4/(4Xh + ^He) is the mean molecular weight of 
the gas in units of the mass of the H atom, mu = 1 .674 lO"^'* g, 
and Hgas is the number density of the gas particles (in cm"^). 

Finally, inserting Eqs. (fTTT ) and ( fT2b into Eq. ( fTOb . we can 
re-write Eq. (|9]) to show the dependence on the dust-to-gas ratio 



3 Agas mn ,„ 

0-„as = Y 



4 ap. 



(13) 



grain 



It is worth noting that in the Elitzur & Ivezic (2001) formal- 
ism the dust-to-gas mass ratio is indeed a free parameter In 
Sect. 14.31 we detail our scheme to specify *!* as a function of 
basic stellar parameters like M and the C/O ratio. 

Moreover, given the dust chemical type (either amorphous 
carbon or silicates) in Elitzur & Ivezic (2001) the quantities Qp 
and 2* are assigned fixed values corresponding to a particular 
choice of the dust condensation temperature (T^ - 800 K) and 
the stellar effective temperature (Tf^g - 2500 K). We have car- 
ried out additional runs of DUSTY (Ivezic et al. 1999) in order 
to evaluate Qp and 2* for various values of T^ and T^s, and for 
different dust compositions. The results are presented in Table[3] 
and are used in our calculations by interpolating in Teff and T^ 
and allowing for different dust mixtures. In the present paper 
we adopt T^ = 1000 K when C/0< 1 and T^ = 1500 K when 
C/0> 1. The grain size is set a = 0.1 fim in all cases. The other 
input parameters, i.e. effective temperature and dust chemical 
type, are just those predicted by our stellar isochrones. 

It should be also mentioned that Eq.(|5]) is valid as long as 
dust condensation is efficient enough to drive an outflow. As 
discussed in Elitzur & Ivezic (2001) this condition corresponds 
to the minimum mass-loss rate: 



M„ 



3 10" 



M^ 



2*0"22^4r^3 



1/2 



(Moyr"') (14) 



where 7^3 is the kinetic temperature in units of 10^ K at the 
inner boundary of the dust shell. For simplicity we assume 
that the kinetic temperature Ty, coincides with the condensation 



temperature T^ of the dust. Interestingly we have verified that 
Mmin given by Eq. ( fT4l l is in good agreement with the results 
from hydrodynamical wind calculations, i.e. Mmin ~ 3 -5 10"^ 
Mq yr"' as predicted by Wachter et al. (2002) for C-stars. 

For M < Mmin, while dust is not a driving mechanism to 
produce a wind, it may be still formed passively in the cir- 
cumstellar ejecta. To handle this case we have simply evaluated 
the expansion velocity setting M = Mmia in the corresponding 
Eq. Q. In this way Vexp is seen to decline smoothly with the 
luminosity to the lowest measured values and, perhaps more 
importantly, by using Vgxp = Vg^piM^in) together with the ac- 
tual M in Eq. Q, we can deal with the very optically-thin cir- 
cumstellar envelopes (T(l//m) ^ 10"^ - 10"-') expected at low 
M. 

One important feature of Eq.(|5]l is that it naturally predicts 
a maximum value for the expansion velocity, which is attained 
at large M as a consequence of the weakening of radiative cou- 



pling in optically thick winds: v„ 



p(tv = 1.3). This 



aspect may bring along significant implications for the inter- 
pretation of the observed data (Sect. lSTl see Fig.|4]l. 

4.3. Dust-to-gas mass ratio 

Under the assumptions of steady-state, spherically symmet- 
ric, radial outflow with complete coupling between dust and 
gas (i.e. both components sharing the same velocity), from the 
mass-conservation law M = An p- pv^^^, the dust-to-gas mass 
ratio ^ - defined by Eq. (fTOl i - can be expressed 



>P: 



Mdu 



M, 



(15) 



gas 



where M^ust and Mgas - M - Mdum denote the mass-loss rates 
in the form of dust and gas, respectively. For any given com- 
bination of M and C/O, Mjust is calculated on the basis of the 
formalism presented by Ferrarotti (2003) and Ferrarotti & Gail 
(2006, and references therein), to whom we refer for many de- 
tails. The general scheme is the following: 



Mdust = ^ Mdust,; 



(16) 



where the summation includes the contribution of various dust 
species, each one being expressed 



A^dust,/ - M Xseed | —. '- I /dust,/ 



(17) 



where Xseed is the total mass fraction of the seectj element in 
the circumstellar envelope, Aseed denotes its atomic weight, 
whereas Adust,; corresponds to the mean molecular weight of 
the dust species. The adopted values are taken from Ferrarotti 
& Gail (2002). The quantity /dust,/ is the condensation degree 
of the dust species, being defined as the fraction of the seed el- 
ement condensed into dust grains. In our calculations the con- 



^ Following Ferrarotti & Gail (2006), with seed or key element we 
denote the least abundant element among those required to form the 
dust species under consideration; e.g. Si is the seed element of the 
silicate compounds in 0-rich stars. 
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Table 3. Characteristic quantities used in the dust models, following the formalism developed by Elitzur & Ivezic (2001). They 
are defined in Sect. 14.21 



dust composition 



Qv PgramCgrcm ) 2i ijim) T^ff Q*{Tss) Qp{T^s) 



100% AlOx 



0.370 



3.2 



0.1 



2000 
2500 
3000 
3500 



0.0189 
0.0344 
0.0599 
0.0958 



0.0087 
0.0105 
0.0128 
0.0157 



100% Silicate 0.779 3.0 



0.1 



2000 
2500 
3000 
3500 



0.0688 
0.106 
0.159 
0.231 



0.0531 
0.0696 
0.0884 
0.110 



1500 

1200 

1000 

800 



100% Silicon Carbide 5.17 



3.2 



0.1 



2000 
2500 
3000 
3500 



0.148 
0.293 
0.527 
0.847 



0.0825 
0.120 
0.176 
0.249 



1500 

1200 

1000 

800 



100% AMC 2.53 



1.9 



0.1 



2000 
2500 
3000 
3500 



0.170 
0.274 
0.425 
0.615 



0.132 
0.180 
0.242 
0.313 



1500 

1200 

1000 

800 



QpiT,) 



1500 0.00811 

1200 0.00897 

1000 0.0107 

800 0.0145 



0.0396 
0.0333 
0.0303 
0.0291 



0.0576 
0.0438 
0.0348 
0.0264 



0.0955 
0.0802 
0.0730 
0.0686 



densation degrees are estimated as a function of the total mass- 
loss rate and C/O ratio by means of the the analytic fits provided 
in Ferrarotti (2003). 

The types of dust species that may condensate in the AGB 
circumstellar envelopes depend critically on the available abun- 
dances of the main dust forming elements, i.e. C, O, Mg, Si, S, 
and Fe. We recall that in the present calculations the chemical 
mixture of metals at the ZAMS is assumed to be scaled solar 
according to the Grevesse & Noels (1993) compilation. While 
the surface abundances of C and O may change in the course 
of the stellar evolution (mainly due to dredge-up episodes and 
HBB), those of Mg, Si, S, and Fe remain practically unaltered. 

Following the results of non-equilibrium dust-formation 
calculations carried out by Ferrarotti & Gail (2006), it is use- 
ful to define two critical carbon abundances, namely Fc,i = 
7o - 2ysi and Yc,i = Fo - ^si + i's, where Y = XjA is the 
abundance in molg '. They define sharp boundaries between 
three different domains of dust compounds, as described in the 
following. 

Adopting the solar O abundance of Grevesse & Noels 
(1993), and under the assumption that the surface abundance 
of O is not affected during the evolution, then Yc,\ and Yc,2 
define two critical C/O ratios, namely (C/0)i = 0.90 and 
(C/0)2 = 0.97. The invariance of O is a common assumption, 
since standard predictions of detailed AGB nucleosynthesis in- 
dicate that the material dredged-up at thermal pulses comprises 
only a small amount of oxygen, i.e. Xq ~ 0.5 - 4% (Boothroyd 
& Sackmann 1988; Izzard et al. 2004). However one should 
consider that this result is highly dependent on the treatment of 
convective boundaries (see e.g. Herwig 2000 for high O inter- 
shell abundances); in addition the surface abundance of O may 
be even depleted as a consequence of eflicient HBB in low- 
metallicity high-mass AGB stars (Ventura et al. 2002). 

Since in our calculations the oxygen variations are negli- 
gible in most cases, we adopt the initial O abundance to de- 
fine three main regimes of dust composition as a function of 
(C/0)i and (C/0)2, which broadly correspond to the three main 



spectral classes of AGB stars (M, S, C). Details are given in 
Sects. l43T1 to l433] 

When applied to our isochrones, the formalism just out- 
lined produces T values typically up to 0.004, 0.002, 0.0001 
for O-rich stars, and up to 0.01, 0.02, 0.02 for C-rich stars, with 
Z - 0.019, 0.008, 0.001, respectively. These are maximum pre- 
dicted values, and they are not meant to be representative of 
AGB stars at varying Z. Detailed population synthesis models 
are needed to explore this issue, which is postponed to future 
analysis. 

Anyhow, from comparing the upper values given above we 
derive some interesting trends. While for O-rich stars we find a 
direct proportionality between ^' and Z, for C-rich stars the re- 
lation becomes more complex and non-mono tonic. In the case 
of O-rich envelopes the relation tfr cc Z is explained by consid- 
ering that the maximum amount of dust that can be formed is 
controlled by the abundances of the seed elements (i.e. Si, Fe), 
which scale with the initial metallicity. 

In the case of C-rich envelopes one key parameter in the 
formation of dust is the excess of carbon with respect to oxy- 
gen, and this quantity does not trace the initial metallicity as it 
is crucially affected by the third dredge-up and HBB during the 
TP-AGB phase. It follows that the large ^' values predicted for 
C-rich stars even at low Z reflect their relatively high C/O ratios 
as predicted by our TP-AGB models at decreasing metallicity. 

It is worth noticing that these theoretical expectations are 
in line with the observational evidence presented in van Loon 
(2000) and Marshall et al. (2004) for the proportionality of the 
dust-to-gas ratio on Z in oxygen-rich stars; and van Loon et al. 
(1999) for the enhanced abundance of free carbon available for 
molecule and dust formation in metal-poor carbon stars. More 
observational data on this issue has been recently provided by 
e.g. Speck et al. 2006, Thompson et al. 2006, Sloan et al. 2006, 
Zijlstra et al. 2006, Lagadec et al. 2007, Matsuura et al. 2007. 
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4.3.1. Stars with C/0 < 0.90 

In this interval there are sufficient oxygen atoms to allow the 
formation of the silicate-type dust. These conditions roughly 
correspond to the domain of M stars. Two main dust contribu- 
tions are considered: 



dt 



+ 



dt 



dt 



(18) 



that correspond to silicate dust and iron dust, respectively. 
These can be explicited 



c/Mci 



dt 
dM„oM 



dt 



Asi 



Ape 



(19) 
(20) 



where Asi (Ape) is the atomic weight of Si (Fe), Asi (Aii-o) is 
the effective molecular weight and /sn (/iro) is the condensation 
degree of the silicate (iron) dust. 

The condensation degree of silicate dust, f^\\ accounts for 
the possible contributions of oUvine-type dust, pyroxene-type 
dust, and quartz-type dust: 



/sil — Jol + /py + Jqu 



(21) 



with clear meaning of the symbols. The eff'ective molecular 
weight of the sihcate dust mixture is expressed 



(./ol Aol + ./py Apy + Jqu Aqu)/jsil 



(22) 



From the analysis of Ferrarotti & Gail (2001), the amount 
of silicate dust in form of olivine, pyroxene and quartz depends 
critically on the Fng/i'si ratio, which is ~ 1.07 according to 
the chemical mixture adopted in this work (Grevesse & Noels 
1993). At this value the non-equilibrium dust condensation cal- 
culations by Ferrarotti & Gail (2001) predict that the dust mix- 
ture is mainly composed by olivine and pyroxene, with a small 
fraction of quartz. Following the results for an M-star model 
with M - 10"^ Mo yr"' and FMg/J'si - 1-07, the condensa- 
tion degrees are assumed to obey the ratios foi/fpy - 4 and 

/ol//qu = 22. 

Finally, /sii and /iio are evaluated with the analytic fits pro- 
vided by Ferrarotti (2003): 



(23) 
(24) 



4.3.2. Stars with 0.90 < C/0 < 0.97 

In this interval a small amount of silicates can be produced 
because of the deficit in O and the dust composition is iron- 
dominated. These conditions may roughly include the class of 
S stars. Two main dust contributions are considered: 




dMA, 



dt 



dM^ilM dMi,oM 

+ 



dt 



dt 



(25) 



that correspond to silicate dust and iron dust, respectively. 
These can be explicited 



dM^ilM .-.^ A sil 

— -— = MXsi--/sii 
dt Asi 

-— = MXpe-— /iro 

dt Ape 



(26) 

(27) 



with the same meaning of the symbols as in Sect. 14.3.11 

Following Ferrarotti (2003) the condensation degrees are 
approximated with 



/sil,S = 0.1 
/iro.S = 0.5 



M_ 



M_6+5 

M-b 
M_6+5 



(28) 
(29) 



We note that for given M the total condensation degree of 
silicate dust is lower than in the case of M stars (because of the 
lower abundance of available O). The relative ratios between 
/ol, /py, and /qu are assumed the same as those for M-stars, 
which is a reasonable assumption (see Ferrarotti & Gail 2002). 

4.3.3. Stars with C/0 > 0.97 

In this interval the dust mixture becomes carbon-dominated. 
These conditions apply to the class of C stars. Two dust contri- 
butions are included: 



t/Mdust liMsicC <iMcai-,C 



dt 



dt 



dt 



(30) 



that correspond to silicon carbide dust and solid carbon dust, 
respectively. These can be explicited 



^Msicc ,-^^ Asic , 
dt Asi 



dMc-. 



car.C 



dt 



= MXcU 



(31) 
(32) 



Xsi, and Xc represent the element abundances (in mass frac- 
tion), which are given by the initial values on the main se- 
quence for Fe and Si, while the current atmospheric value is 
adopted for C. The quantities Asi and Asic are the atomic weight 
of Si, and the effective molecular weight of the silicate dust 
mixture, respectively. 

Finally, /sic, and /car denote the degree of condensation 
of the elements Si and C in the corresponding dust species. 
Following Ferrarotti (2003) they are evaluated with 



/s.c = 0.5 

/car = 0.5 



M. 



M_6+5 

M_6 /C 
M-6 + 5 \0 



(33) 
(34) 



4.4. Radiative transfer calculations 

Bolometric corrections are provided from two independent RT 
calculations, which are updates and extensions of those origi- 
nally described by B98 and G06. In both cases the models have 
been calculated with 1 -dimensional codes that solve the RT 
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equation and the thermal balance equation in a self-consistent 
way (see Granato & Danese 1994; Groenewegen 1993, and 
Groenewegen 1995). 

Although they differ in several minor details which are re- 
lated to the structure of the original codes - as described in 
Sects. 14.4.11 and 14.4.21 - B98 and G06 calculations produce 
very similar results, which proves the robustness of both RT 
models. For instance, whereas B98 tabulate the complete out- 
put F,i spectra with a regular spacing in log A (here set to be 
A log A = 0.04 dex), G06 code provides directly the integrated 
F^.j for a given set of input pass-bands. In fact, we have es- 
timated that the differences in ABC^.i^ coming from these two 
different approaches, for the broad-band filters we are using, 
are of the order of just ~ 0.05 mag forT(lyum) = 1. Considering 
the overall uncertainties, these differences are small enough to 
be safely ignored. 

Much more important, instead, are the differences associ- 
ated with the assumed dust compositions, and the coverage of 
stellar parameters, as detailed below. 



^ 




log K„ (^) 



4.4.1 . Extended G06 calculations 

The dust models that are incorporated in the isochrones are an 
extension of the work described in G06. The models have been 
calculated for (arbitrary) values of L = SOOOL©, li = 8.5 kpc, 
Vexp = lOkms ', T = 0.005. In G06, only few effective tem- 
peratures were considered. Here we calculated models for spec- 
ti-al types M0,1,. . .,10 (reff = 3850, 3750, 3650, 3550, 3490, 
3397, 3297, 3129, 2890, 2667, 2500 K) fromFluks et al. (1994) 
for O-rich stars, and T^g = 3400, 3200 , 3000, 2800, 2650 K 
from Loidl et al. (2001) for C-rich stars (with a C/O ratio of 
1.1). The O- and C- stellar input spectra are for solar metallici- 
ties. 

Several types of dust are considered (see G06 for a 
discussion). For O-rich stars, (1) 100% Aluminium Oxide 
(AlOx; amorphous porous AI2O3), with optical constants from 
Begemann et al. (1997), (2) a combination of 60% AlOx and 
40% Silicate (optical constants from David & Pegourie 1995), 
(3) 100% silicate. A condensation temperature of Tc = 1500 K 
was assumed in the first two cases, 1000 K for the pure silicate 
dust. 

For C-stars, we consider either a combination of 85% 
Amorphous Cai'bon (AMC) and 15% Silicon Carbide (SiC) 
with optical constants from, respectively. Rouleau & Martin 
(1991) for the ACl species and a-SiC from Pegourie (1988), 
and assuming a T^ of 1200 K, or 100% AMC with a T^ of 
1000 K. 

The results of such computations consist in the apparent 
stellar magnitudes for several values of M (and hence t) rang- 
ing from lO"' to 4 X 10"^ Mq yr ', and for each spectral type 
and filter under consideration. Over this range of mass loss t 
increases up to ~ 3 for AMC (at 1 1.33 yum) and ~45 for AlOx 
(at 11.75 fim). A total of 90 filters were considered in these 
specific computations. These magnitude tables were converted 



Fig. 3. The absorption/emission properties of the dust mix- 
tures considered in this work, in the form of ABC values and 
limited to the case of ri^m = 1 and input stellar spectra of 
Tss = 3000 K. Top panel: For O-rich stars, the dust compo- 
sitions are 100% AlOx (crosses), 100% Silicates (triangles), 
60% AlOx + 40% SiHcate (squai-es) following the G06 ap- 
proach, and Silicates (continuous line) following B98. Bottom 
panel: For C-rich stars, the dust compositions are 100% AMC 
(crosses) and 85% AMC +15% SiC (squares) following the 
G06 approach, and Graphites (continuous line) following B98. 



into ABC(t) values that directly give the difference between the 
magnitudes of non-dusty (r = 0) and dusty staro 

4.4.2. B98 calculations 

In B98 for each dust type (either O-rich or C-rich) the refer- 
ence grid of AGB dusty envelopes consists of models calcu- 
lated for L = 10000, Tetf = 3000/: while varying the opti- 
cal depth T(lyum) from 5 10""* to 200 in 50 equal logarithmic 
steps. The condensation temperature is assumed T^ - 1000 K 
for O-rich envelopes, and T^ - 1500 K for C-rich envelopes. 
Additional tests have shown that, for the adopted dust mixtures, 
the resulting spectral shapes are little affected by changes in L 
and Teff within the typical ranges for AGB stars, for fixed val- 
ues of r(l//m). 

The radiative transfer of the photospheric stellar spectrum 
through the dust is computed with the Granato & Danese 
(1994) code (originally built for AGN tori), assuming spherical 
symmetry and the radial dust density distribution decreasing as 
r"^. The adopted dust grain model is that of Rowan-Robinson 



The use of ABC(t) ensures that any residual difference in the 
zero-points between the G06 magnitudes, and those used in our new 
isochrones, will be erased. 
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et al. (1986)13, which considers six different kinds of grains: 
AMC of radius a - O.ljum and graphites of 0.03 and 0.01/im 
for C-rich mixtures, and amorphous silicate of a = 0. 1/im and 
silicates of 0.03 and 0.01/vm for O-rich mixtures. 

Although we updated in several aspects the B98 model, in 
particular in the dust grain properties, the detailed exploration 
with Spitzer data of these updates is still underway and will be 
presented in a forthcoming paper. Therefore here we adopt the 
already well established B98 model in its original form. The 
dust composition we consider is carbon for C-type mass-losing 
stars, and silicates for M-type ones. 

The relation between T(ljum) and the mass loss is given by 
Eq. dUl with Ad = 1.5 X 10" for silicates and Ad = 6.6 x 10" 
for the carbonaceous mix (with Vexp in km s" 1, L in solar units). 
The quantity Ad is slightly affected by the spectrum of the il- 
luminating star. However, B98 have checked that the errors in 
T(l//m) estimated from Eq. (|4]i with the above values of Ad are 
kept within 5% for spectra with Tetf in the range 2500-4000 K. 
This is a fair approximation because models with T(l/im) dif- 
fering by less than 5-10% produce almost identical spectra. 




1000 2000 

P (days) 



1000 2000 

P (days) 



4.4.3. Overview of RT models 

Figure|3]illustrates the behaviour of ABC as a function of wave- 
length, for the dust compositions considered in this work, and 
for the specific case of T(lyum) - 1 and T^s - 3000 K. It can 
be noticed that all curves agree very well in the case of O-rich 
mixtures, and differ significantly only for the C-rich mixtures at 
A > 10/Ym. These differences become higher for higher values 
of T, and can be ascribed to the use of different dust mixtures 
and spectral libraries. 

In summary, we consider several possible choices of dust 
composition, namely: four for O-rich envelopes and three for 
C-rich envelopes. Alternatively, we can ignore dust setting r = 
0. All these cases are allowed in our isochrone calculations. 

5. A preliminary comparison with observations 

GiOO isochrones have been extensively used in the literature for 
the study of resolved stellar populations. Most of their use has 
regarded the optical photometry of stars in the most populated 
evolutionary phases like the main sequence, sub-giant branch, 
RGB, red clump and early- AGB. Now the addition of the de- 
tailed TP- AGB evolution together with the infrared photometry 
open new perspectives for the application of these new stellar 
models, that we briefly illustrate in the following via compar- 
isons with Galactic Disk and Magellanic Cloud data. 

Most of the comparisons are intended to illustrate the lo- 
cus of our new isochrones in different diagrams, and do not 
consider crucial aspects of the host galaxies such as their 
age-metallicity relations and star formation histories. A de- 
tailed comparison between present models and the observed 



' The dust model by Rowan-Robinson et al. (1986) is based on ob- 
servations of circumstellar dust shells around late-type C- and M-type 
stars, that appear to require absorption efficiencies decreasing as A^' 
in the far-IR. This wavelength dependence is ascribed to grains with 
a disordered, damaged or amorphous structure. See the original paper 
for more details. 



Fig. 4. Expansion velocities of circumstellar flows against pul- 
sation periods of variable AGB stars, for both O-rich (left pan- 
els) and C-rich objects (right panels). Upper panels: Galactic 
observations (empty circles) are taken from the compilations 
assembled by Groenewegen et al. (1998) for O-rich stars and 
Groenewegen et al. (2002) for C-rich stars. Observed data for 
LMC variables (Marshall et al. 2004) is also plotted (empty 
squares). The fit relation proposed by Vassiliadis & Wood 
(1993) is shown by the solid line. Bottom panels: Predicted ex- 
pansion velocities for a few selected isochrones with Z - 0.019 
and Z - 0.008. Ages are indicated in the figure. See the text for 
more details. 

star counts in several parts of these diagrams, is left to future 
papers. 

5.1. Expansion velocities of variable AGB stars 

As first check of our models we compare the predicted expan- 
sion velocities (Sect. I4.2l i as a function of pulsation periods 
with observations of variable AGB stars in the Galactic Disk 
and Magellanic Clouds. The results, displayed in Fig. |4] are 
consistent with the empirical estimates. 

A relevant remark should be made in respect to the O-rich 
stars. While for P < 500 - 600 days a positive correlation be- 
tween periods and expansion velocities seems to exist in the 
observed data, this trend does not proceed at longer periods, 
where a considerable scatter of Vexp is detected. It should be no- 
ticed that while reproducing the mean observed trend at shorter 
periods, the widely-used relation of Vassiliadis & Wood (1993; 
solid line) fails to describe the data at longer periods, where 
Vexp is assumed to be constant (15 km s"'). 

In the bottom panels of Fig. |4] we show the results of our 
calculations for a few selected isochrones with Z = 0.019 and 
Z - 0.008. Regardless of age all isochrones describe a similar 
non-monotonic evolution on the Vgxp -P diagram. At increasing 
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periods, first the expansion velocity increases, reaches a max- 
imum and eventually declines. This bell-shaped behaviour is 
particularly well developed for the youngest isochrones (with 
ages of ~ 1.6 10** yr and ~ 2 10^ yr), which correspond to mas- 
sive AGB stars (with initial masses of 4.5 Mq and 4.1 Mq) and 
it appears to explain the existence of variables with high peri- 
ods and relatively low expansion velocities. 

To this respect it is also interesting to compare the two 
isochrones with log(age/yr) = 8.3) and different metallici- 
ties (Z = 0.019 and Z = 0.008). The one with Z = 0.008 
reaches, on average, lower Vexp- Moreover, the few OH/IR stars 
in the LMC for which Vexp has been measured (empty squares; 
Marshall et al. 2004) fall just in correspondence of the declin- 
ing branches of the v'exp curves. This fact suggests that the low 
expansion velocities of LMC OH/IR stars compared to those of 
Galactic OH/IR stars with the same P, may not reflect a pure 
effect of metallicity (i.e. a linear Z-dependence of the dust-to- 
gas ratio *P such that the lower Z the lower Vexp), as commonly 
invoked (see the discussion in Marshall et al. 2004); rather it 
may be the result of a more complex dependence of v'exp as a 
function of stellar and dust parameters (L, M, M, *F). A deep 
investigation of this issue will be performed in future work. 

5.2. The LMC and SMC near-infrared photometry 

TP-AGB stars become striking objects in galactic surveys con- 
ducted at pass-bands redder than /. The best available data for 
the AGB population in galaxies is likely the one obtained from 
the extensive near- and mid-infrared surveys of the Magellanic 
Clouds. 

In Fig.|5]we plot the 2MASS Ks vs. J-K^, diagram of the in- 
ner LMC and SMC. The data are taken directly from the SAGE 
and S^MC catalogs (Blum et al. 2006, and Bolatto et al. 2007, 
respectively), and hence correspond to objects that have coun- 
terparts in the Spitzer pass-bands. 

We start discussing the LMC case depicted in the left panel. 
The circular area of tt sqrdeg centered on the LMC bar at 
(22000 = 5''23™5,52ooo = -69°45', is large enough to include 
thousands of AGB stars brighter than the TRGB, and small 
enough to somewhat limit the presence of features that are not 
due to the LMC. These latter features are however still evident 
in the plot, and consist of (1) the vertical fingers in the top left 
part of the CMD, at K^ < II mag and J-K^ < 1 .0 mag, caused 
by foreground stars in the Milky Way disk and halo, and (2) the 
plume of objects at the bottom right part of the CMD, roughly 
at ^s > 12 mag vs. J-Kf^ > 1.6 mag, which largely consists 
of background galaxies (see Nikolaev & Weinberg 2000, his 
Region L). Outside of these two CMD regions, the bulk of ob- 
jects indeed consists of stars located in the LMC. 

On this diagram, we over-plot our isochrones for the ini- 
tial metallicity Z - 0.008 (which are adequate to describe the 
youngest LMC populations) and at four selected ages: 

- At log(f/yr) - 7.3, AGB stars are not present and the most 
prominent near-IR stars are red super-giants burning He in 
their cores. These stars draw a diagonal line which coin- 
cides with sequence H of Nikolaev & Weinberg (2000; their 
Fig. 3). 



- At log(f/yr) = 8.2, the TP-AGB phase is well developed 
and populated by luminous O-rich stars undergoing HBB. 
These stars draw a second diagonal sequence in the CMD, 
corresponding to sequence G of Nikolaev & Weinberg 
(2000) and joining their sequence F at fainter magnitudes 
which is populated by early-AGB stars. These features are 
also well described by our new isochrones. 

- At log(f/yr) - 9.1, HBB does not occur and the TP- 
AGB phase is dominated by a well-populated red tail of 
C-rich stars (sequence J in Nikolaev & Weinberg 2000), 
which departs from the O-rich sequence at K^^ = 10.3 mag, 
J-Ks - 1.2 mag and reaches J-Kf^ as red as 2.0 mag. 
This feature is overall well described by the log(f/yr) = 9.1 
isochrones without dust, apart from the fact that in this case 
the red tail just reaches J-K^ = 1.6 mag, which is the 
colour of all C-star spectra cooler than ~ 3000 K in the 
Loidletal. (2001) library. 

- At log(f/yr) = 10.0, the third dredge-up does not occur, 
and the only TP-AGB stars are O-rich. They correspond to 
sequence F in Nikolaev & Weinberg (2000), located right 
above the TRGB. 

The above items mainly describe the behaviour of the 
isochrones at the regime of weak mass-loss, in which colours 
and magnitudes are affected by circumstellar dust at the level of 
just $ 0.2 mag. What happens in the regime of high mass-loss 
rates (say \M\ > 10"^ Mq yr"') is also well evident in the figure, 
when we look at the isochrones with dust obscuration: the up- 
per fraction of their TP-AGB phase departs to very red colours 
at increasing /STj-band magnitudes, reaching J-K^; > 4 mag. 
The location and slope of this extended red tail (sequence K 
in Nikolaev & Weinberg 2000) is also well described by our 
isochrones, and especially by the log(f/yr) = 9.1 one which 
runs over the most densely populated part of this extended red 
tail. 

The left panel of Fig. |6] details the region of the K^ vs. 
J-K^ diagram that includes the bulk of the TP-AGB stars. In 
addition to the 2MASS data for the inner LMC, we plot the 
Groenewegen (2004) database of LMC stars spectroscopicaUy 
classified as being of types M and C. The latter are LPVs se- 
lected from the OGLE database for the entire LMC. Although 
they are not distributed evenly over the CMD, they clearly il- 
lustrate that the bulk of M- (dark-green) and C-type (orange) 
stars are well separated in the CMD; in particular they show 
the striking red tail of C stars extending from J-Kf^ ~ 1.3 mag 
up to y-ZTs ~ 2.0 mag, and the confinement of M stars at the 
bluest colours (J-Kf^ ^1.3 mag). 

Then, the figure evidences how the new isochrones repro- 
duce the general appearance of the data, and in particular the 
red tail of C stars which is well described by the log(f/yr) = 9.1 
isochrone. It is also remarkable how well the isochrones pre- 
dict the neat separation between O-rich (M-type) and C-rich 
(C-type) stars in the diagram. In the log(f/yr) = 9.1 isochrone, 
the red tail reaches J-Kf^ ~ 1.6 even in absence of circum- 
stellar dust. Dust obscuration in stars with \M\ 5 10"^ Mq yr"' 
provides an additional A(J-Ks) ~ 0.2 - 0.3 mag that improves 
the comparison at the reddest part of the red tail, at J-K^ be- 
tween 1.6 and 2.0 mag. As thoroughly discussed in Marigo et 
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Fig. 5. The 2MASS Kg vs. /-/Ts diagram for the inner LMC (left panel) and for the SMC (right). Over-imposed are isochrones 
with Z - 0.008 shifted by a distance modulus of 18.5, and with Z - 0.004 shifted by 19.0, for the LMC and SMC cases, 
respectively. The isochrones are drawn in different colours according to their evolutionary status: pre-TP-AGB stars (green), 
and TP-AGB stars O-rich (blue) and C-rich (red). The selected ages are: log(f/yr) = 7.3 (no AGB present), 8.2 (well de- 
veloped O-rich TP-AGB), 9.1 (well-developed C-rich TP-AGB), and 10.0 (short O-rich TP-AGB), for 3 different choices of 
dust-obscuration: the no-dust case (continuous lines), cf. Bressan et al. (1998; dot-dashed lines) and c.f. Groenewegen (2006; the 
0.85 AMC+0. 15 SiC mixture for C-stars, and 0.6 AlOx+0.4 Silicate for O-rich stars; dashed lines). In order to limit confusion in 
this plot, the isochrones are drawn up to their point of maximum mass-loss rate on the TP-AGB, i.e. ignoring their final post- AGB 
stages in which the isochrones cross back to the blue part of the diagram. Moreover, a magenta dot mark points where the mass 
loss increases above 10"^ Mq yr"'; above this indicative limit the stellar lifetimes are necessarily short, and few stars are expected 
to occupy isochrone sections to the red of these points. 



al. (2003a), the main physical effect driving the appearance of 
the red tail of C stars are the cool Teg caused by changes in 
molecular opacities as the third dredge-up events increase the 
photospheric C/O ratio. 

The right panel of Fig. |6] over-plots isochrones from GiOO 
and Cordier et al. (2007) over the same data. These isochrone 
sets contain the complete TP-AGB phase but computed in a 
rather crude way, i.e. without considering the third dredge-up 
events and HBB nucleosynthesis, and using Teff (L, M, Mcoie) 
relations derived from scaled-solar mixtures. These isochrones 
fail in reproducing the basic features indicated by 2MASS and 
spectroscopic data, although in different ways: 

- The GiOO isochrones describe well the locus of O-rich AGB 
stars, as expected for TP-AGB models that do not consider 
the third dredge up process. These isochrones are confined 
to J-Kfi $ 1.3 mag, and do not cross the red tail of C stars 
between J-Kf^ ~ 1.3 and 2.0 mag. 

- The Cordier et al. (2007) isochrones, instead, do produce 
a sort of red tail, reaching J-K^^ ~ 2.1 mag, but they do 
this while still on an O-rich phase, and for all ages younger 
than about 1 Gyr. This is completely at odds with obser- 
vations, which show that only C-rich stars (with very few 
exceptions) populate the red tail at J-K^ > 1.3 mag. 



The reasons why Cordier et al. (2007) O-rich TP-AGB models 
do reach so red J-K^ colours - so as to mimic C-type stars, 
though they are not - likely reside in both the extremely low 
Teff - as cool as 2400 K at the TP-AGB end - they reach using 
the Wagenhuber (1996) formula, and on the use of a colour- Teff 
relation derived from C-type stars (from Bergeat et al. 2001). 
These choices are highly questionable. First, independent inte- 
grations of model atmospheres (e.g. Marigo 2002) suggest that 
O-rich stars can hardly reach Tgff values lower than ^ 3000 K, 
even at solar metallicitieo- Second, it has long been known 
that (J -K)-T f^ff relations for C-stars cannot be applied to O- 
rich stars, and vice-versa, since these two kinds of stars have 
very different near-infrared spectra. 

From the failure of both the GOO and Cordier et al. (2007) 
isochrone sets in describing a red tail made of C stars, it follows 



'" Note that the bulk of M giants in the Magellanic Clouds are of 
spectral types earlier than M4, and in general much hotter than the typ- 
ical M giants in the Milky Way. A handful of cool M giants and super- 
giants of spectral type later than M6 is also known in the Magellanic 
Clouds (e.g. Groenewegen & Blommaert 1998; van Loon et al. 1998, 
2005). In order to explain them, and also the more common giants of 
spectral type later than M8 observed in the solar neighbourhood (with 
Teft- as small as 2400 K, see Fluks et al. 1994), our atmosphere models 
require super-solar metallicities. 
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Fig. 6. Left panel: The 2MASS Kf^ vs. J-Ks diagram for the LMC, now detailing the CMD region that corresponds to the bulk 
of TP-AGB stars, and limiting the isochrones to the ages of log(f/yr) = 8.2, 9.1, and 10.0. Additional LMC data are also plotted: 
The small circles are LPVs in the LMC that have been spectroscopically classified as being of type M (dark-green) or C (orange), 
according to the data compilation from Groenewegen (2004). The meaning of the different lines along the isochrones are the 
same as in Fig.|5] The right panel shows the same data as compared to GiOO (blue lines) and Cordier et al. (2007; magenta lines). 
Notice that these latter sets of isochrones do not contain a C-star phase. 



that such models are not suitable for the interpretation of near- 
infrared photometry in the age interval in which C stars form. 
For LMC metallicities, this interval is quite ample, going from 
a few 10^ yr up to about 5 Gyr (see also Frogel et al. 1990). 

Finally, the 2MASS ^s vs. J-K^ diagram for the area 
within the SMC covered by the S^MC survey (right panel of 
Fig. |5]l can be described in quite a similar way to the LMC 
one. The most obvious difference between the two panels is 
in the ~ 0.5 mag shift of all features coming from the LMC- 
SMC difference in distance, that we also took in consideration 
when plotting the isochrones. Also, due to the particular se- 
lection of data the SMC plot presents less contamination from 
background galaxies than the LMC one. 

5.3. The LMC and SMC mid-infrared piiotometry 

Figure s|7] and |8] show a series of CMDs involving the near- and 
mid-IR photometry of the Magellanic Clouds, comparing the 
new isochrones with 2MASS+SAGE data for the inner LMC 
(Blum et al. 2006) and with 2MASS-hS^MC data for the SMC 
(Bolatto et al. 2007O In all these plots, the effect of circum- 
stellar dust becomes crucial in describing the observed pho- 
tometry. The features present in these CMDs (except for the 
IRAC-hMIPS [24] vs. [8.0] - [24] diagram, for which the data 
is poorer) are essentially the same as those present in the K^ vs. 
J-K^ diagram. In particular, all of them show contamination 



" The S^MC data, originally in flux units, have been converted to 
magnitudes using the zero-points from Cohen et al. (2003) for JHK^ 
and from Bolatto et al. (2007) for Spitzer pass-bands. 



from foreground Milky Way stars in the form of a narrow ver- 
tical finger confined to the left part of the CMDs, and contam- 
ination from background galaxies at their bottom-right comer. 
The gray dots correspond to objects not detected by 2MASS, 
most of which are expected to correspond to background galax- 
ies and young pre-main sequence stars, which are not of interest 
in this paper. 

The LMC population in these plots again shows 0-rich 
sequences (core He burning red super-giants plus 0-rich TP- 
AGBs) that progressively merge into an almost-vertical blue 
sequence when colours involving redder pass-bands are con- 
sidered; plus an inclined C-rich sequence that can be identified 
as the main body of the J-Kf^ < 2.0 mag red tail. From the 
extremity of this more populated red tail a very extended and 
less populated plume departs, that corresponds to the stars with 
high mass-loss rates. Our isochrones reproduce well the posi- 
tion of these plumes in the several CMDs, and moreover predict 
that they are populated by both O- and C-rich stars. 

5.4. Integrated broad-band colours 

Broad band colours of Magellanic Cloud clusters represent the 
basic observable against which any population synthesis model 
should be compared for testing its predicting performance. 

Figure |9] displays the results for two integrated colours, 
V-K and J-K, that are expected to be significantly affected by 
the existence of TP-AGB stars. By looking at Fig. |9] (top pan- 
els) we actually see that the observed data for both V-K and 
J-K colours show a rising trend starting from f ~ 10^ yr to- 
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Fig. 7. A comparison between the new isochrones and the near- to mid-infrared data for the LMC. The isochrones are the same 
as plotted in the left panel of Fig. |5] The data (dots) come from the SAGE epoch 1 point source catalog (Blum et al. 2006) and 
combine IRAC+MIPS with 2MASS data inside a circle of 1 degree radius in the inner LMC. The panels show the 2MASS-I-IRAC 
[3.6] vs. y- [3.6], 2MASSH-IRAC [8.0] vs. 7- [8.0], IRAC [8.0] vs. [3.6] - [8.0], and IRACh-MIPS [24] vs. [8.0] - [24] diagrams. 
The gray dots correspond to objects not detected by 2MASS, most of which are expected to correspond to background galaxies 
and young pre-main sequence stars, which are not discussed in this paper. 



wards older ages, which is currently well explained by our cur- 
rent single-burst stellar populations (SSP) models (solid line) 
as due to the development of the TP-AGB phase. The reader 
is referred to Appendix lAl for a discussion on the behaviour of 
V - K ?A ages t ~ 10^ yr and its relation to the evolution of the 
most massive AGE stars. Integrated magnitudes are calculated 
accounting for the luminosity variations (i.e. the long-lived lu- 



minosity dips in low-mass stars) driven by thermal pulses, and 
the effect on the emitted spectrum by circumstellar dust. 

For comparison we also show the results of other two 
widely-used SSP models, namely Chariot & Bruzual (2003; 
dot-dashed line) and Maraston (2005; dashed line). A few 
points should be noticed in respect to the age interval t ^ 
10^ yr, when AGB stars are present. While Chariot & Bruzual 
(2003) colours tend to be bluer than the observed ones be- 
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Fig. 8. As Fig.|2]but for the SMC. The data (dots) come from the S^MC catalog (Bolatto et al. 2007). The isochrones are the 
same as in the right panel of Fig. |5] 



cause of a likely underestimation of the TP-AGB phase (see 
Bruzual 2007 for a recent revision of this aspect in their mod- 
els), Maraston (2005) tends to form a sort of upper envelope 
to the data, especially in V - ^. Our results are intermediate 
between the two. 

The bottom panels of Fig. |9] illustrate the same colour evo- 
lution, by means of simulated clusters. The cluster simulations 
are performed by randomly adding stars to the isochrones, with 
the probability of a given initial stellar mass being given by 
the Chabrier (2001) log-normal initial mass function. The ad- 
dition of stars to a simulated cluster is stopped when its initial 
total mass reaches a given value Mj. The entire cluster sam- 



ple is assumed to follow a M;^:'^ distribution, which is derived 
from the luminosity function of LMC clusters (Elson & Fall 
1985; see also Girardi et al. 1995, sect. 8.2). The cluster metal- 
licities follow the LMC age-metallicity relation from Pagel & 
Tautvaisiene's (1998) bursting model, which has been shown to 
reproduce the observed trend. The distribution of cluster ages is 
assumed to be constant; although this is not a realistic assump- 
tion (see Girardi et al. 1995), it allows us a good illustration 
of the colour evolution over the complete age range. A poste- 
riori, we impose a limiting magnitude of MJ^™ - -6.5 mag to 
the simulated sample, which roughly corresponds to the one of 
the observed samples in the top panels of Fig. |9l The results 



20 



Marigo et al.: Optical to far-IR isochrones 



4 
3 
2 
1 

4 r 
3 r 
2 ^ 
1 - 






L^g^A^ ' 



8 9 

log(t/yr) 



10 



I 



1.5 
1 

0.5 


1.5 
1 

0.5 





1 •• ,^'~ •" •* 




1 1 1 1 1 




Fig. 9. Broad-band V - K and J - K colours of LMC clusters as a function of cluster age. Top panels: Observed photometric 
data taken from various compilations: Persson et al (1983, filled circles); Kyeong et al. (2003, empty squares); and the V-band 
photometry by Goudfrooij et al. (2006) adopting an aperture radius of 50 arcsec, combined to the JHK^ photometry from Pessev 
et al. (2006; starred symbols). Clusters are assigned ages by adopting the 5 -parameter-age calibration by Girardi et al. (1995). 
Predicted SSP colours from various authors are superimposed for comparison: this work (solid line, Z = 0.008), Maraston (2005; 
dashed line, Z = Zq/2), Bruzual & Chariot (2003; dot-dashed line, Z = 0.008). We note that the spike in the colour evolution 
at ~ 10^ yr is real and is explained in detail by Girardi & Bertelli (1998). Bottom panels: Simulations of stellar clusters based 
on the new isochrones. They roughly follow the distribution of typical sizes (or initial total masses) and metallicities of LMC 
clusters. The age distribution is assumed to be constant. The simulated sample has a limiting magnitude of M^" = -6.5, which 
roughly corresponds to the one of the observed samples in the top panel. 



present the intrinsic colour dispersion expected for the clus- 
ter data due to their small numbers of evolved stars (stochastic 
fluctuations), and in this aspect they look clearly more similar 
to the observations than the continuous narrow SSP lines de- 
picted in the top panel. The main discrepancy between models 
and simulations is at ages younger than 10^ yr, where a bunch 
of very blue clusters is predicted but absent from the observed 
samplq^. The blue clusters are essentially those for which no 
red super-giant was predicted, and derive naturally from the 
fact that the red super-giant phase is short-lived. In order to 
recover the observed colours of the youngest clusters the red 
super-giant phase has to be significantly longer and cooler than 
in the Bertelli et al. (1994) isochrones. On one hand, a longer 
duration could be attained, for instance, by invoking a larger 
efficiency of overshoot, but this may in turn inhibit the devel- 
opment of the blue loops of massive stars (see the discussion 
in Alongi et al. 1991). On the other hand, the effective temper- 
atures of red super-giants depend very much on the treatment 
of the density inversion that is predicted to take place in the 
outermost atmospheric layers. A deeper analysis of this issue 
is beyond the scope of the paper and it is postponed to future 
work. 

Another questionable point involves the ages between 
10*^ yr and ~ 4 x 10** yr, after the AGB phase develops: in 



Isochrones in this age range were taken from Bertelli et al. (1994). 



this age interval, our models become on the mean redder than 
the data. As mentioned in the Appendix, removing this feature 
from the models poses a theoretical difficulty as it would as 
well prevent the significant enrichment of helium and nitrogen 
in the most massive AGB stars, which is instead required to ex- 
plain the observed abundances in planetary nebulae of type I. 
Finally, at older ages the adoption of an age-metallicity relation 
allows us to reproduce the colours of the oldest LMC clusters, 
which have ages and metallicities typical of old Galactic glob- 
ular clusters. 

In short, these isochrones offer a valid alternative for the 
modelling of the integrated light of star clusters and galaxies, 
since they pass the LMC cluster test performing similarly well 
as the widely used Chariot & Bruzual (2003) and Maraston 
(2005) SSP models, when we look at the mean colour evolution 
(top panels of Fig.|9]l. Moreover, the isochrones we provide do 
allow the simulation of resolved stellar populations, which is 
essential for their further testing - as illustrated by the simula- 
tions in the bottom panels of Fig.|9] The isochrones are based 
on internally consistent sets of tracks, which naturally obey the 
fuel consumption theorem, and provide other quantities neces- 
sary for a consistent modelling of galaxies - namely the chem- 
ical yields and remnant masses (see Marigo & Girardi 2001 
for a discussion on this latter point). Having all these compo- 
nents together allows us the simultaneous comparison of the 
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model results with many different observables, indicating hid- 
den problems (like the young blue clusters predicted in Fig.|9]l 
and hopefully pointing the way to future improvements in the 
models. 



6. Data tables and final remarks 

6.1. Data tables 

The data presented in this paper are provided in two different 
ways: 

(1) As a repository of static tables down- 
loadable from the VizieR Catalogue Service at CDS 
( http : //vizier . u-strasbg . fr ) and from the URL 
|http: //stev. oapd. inaf ■ it/dustyAGBS7[ This database 
includes not only the isochrones, but also bolometric cor- 
rections, extinction coefficients, dust attenuation curves, 
etc. 

(2) Via an interactive web interface at 
|http : //stev . oapd . inaf. it/cmd that allows users to 
build isochrones (or sequences of them) for different values of 
age and metallicity, photometric systems, and dust properties. 
The same web form provides access to a series of other ser- 
vices, like the computation of luminosity functions, integrated 
magnitudes, and simulated star clusters. 

All information about these data are provided in the form 
of either ReadMe files in case (1), or HTML help pages in case 
(2). Most of the information available will be of interest only 
to the people who actually access the data, and therefore will 
not be repeated here. Suffice it to recall that isochrones are pro- 
vided for any metal content in the range 0.0001 < Z < 0.03 
(corresponding to populations with -2.31 < [M/H] < -hO.22, 
for scaled-solar initial compositions), and ages in the range 
7.8 < log(f/yr) < 10.2 (i.e. from 0.063 to 15.8 Gyr). The initial 
mass ranges from 0.15 to 7 Mq. 

Moreover, the |http: //stev .oapd. inaf .Tt/ cmd URL 
allows the mass range to be extended from 7 to 100 Mq, us- 
ing the tracks from Bertelli et al. (1994) isochrones. This al- 
lows the age range to be extended down to about 4 Myr for 
complete isochrones - or even down to f = in the case one 
is interested in obtaining an extended zero-age main sequence. 
The only problem with this is that the tracks for massive stars 
are not available for Zq < Z < 0.03, and in this case have to be 
taken from the Z = 0.02 set of Bertelli et al. (1994). 

The basic information provided along isochrones consists 
on the luminosity, effective temperature, surface gravity, and 
absolute magnitudes in the photometric system under consid- 
eration. For the TP-AGB phase, this information is comple- 
mented with the core mass, surface C/O ratio, pulsation mode 
and period, and mass-loss rate. The TP-AGB section of an 
isochrone can be easily identified by the information about the 
core mass, whereas the C-star phase is identified by C/O > 1. 
These latter variables are also useful to simulate, a posteriori, 
the luminosity and Tetf variations caused by pulse cycles on the 
TP-AGB, as we have done in a series of papers (e.g. Marigo et 
al. 2003; Cioni et al. 2006; Girardi & Marigo 2007b). 

Finally, we remark that the 

|http : //stev . oapd . inaf . it/ cmd web interface is go- 



ing to be continuously expanded, either with new extensions 
to the tracks (e.g. with the planetary nebulae and white dwarf 
sequences), new photometric systems, new spectra from 
which to derive bolometric corrections, and different dust 
compositions. Interstellar extinction will be included in the 
isochrones in a consistent way (Girardi et al., 2007, in prep.). 
In any case, the original database, as here described, will 
remain available as well. 

6.2. Final remarks 

Users of the new isochrones should be aware that: 

- The isochrones are expected to perform reasonably well 
for AGB populations in the Milky Way disk and in the 
Magellanic Clouds, against which data they have already 
been compared or calibrated. The quality of the results at 
very low metallicities ([Fe/H]^ -1.5) is yet to be verified. 

- Present-day model atmospheres for cool stars are still far 
from satisfactory. Molecular line lists are acknowledged to 
be incomplete, whereas for LPVs dynamical model atmo- 
spheres would probably provide more realistic spectral en- 
ergy distributions than the static models we have adopted 
for both M- and C-type stars. The template spectra used in 
or models do not cover well enough the expected range of 
metallicities, effective temperatures, and surface gravities 
of AGB populations in resolved galaxies. 

- Dust properties may differ from the compositions assumed 
in this work, and they likely change along the TP-AGB 
evolution as a function of parameters other than those here 
considered. For a given dust composition, RT models have 
been calculated for a large range of optical depths, even 
if particular combinations of dust and mass-loss rate are 
unlikely. For example, pure AlOx is preferentially seen in 
stars with low mass-loss rates only, and the combination 
of AlOx and silicates is most appropriate for intermedi- 
ate mass-loss rates. For the Galaxy, the mass-loss rate be- 
low which the pure AlOx models should probably be used, 
and above which pure silicates are likely to be most appro- 
priate are, respectively, about 10"^ Mq yr~' and 10"^ Mq 
yr"' . At different metallicities these ranges may be different 
and therefore the full range of models has been calculated 
for all compositions here considered. An extended list of 
other caveats associated with dusty circumstellar envelopes 
is presented in G06. 

- Mass-loss prescriptions from evolved stars are still very un- 
certain, and probably none of the presently available for- 
malisms can be applied with confidence. For O-rich AGB 
stars, even the basic physical process driving mass loss is 
being debated (see Woitke 2006). 

- In constructing our theoretical isochrones, we have made 
all possible effort to ensure internal consistency, but this 
has not always been possible. For instance, the dust con- 
densation prescriptions here adopted are likely inconsistent 
with those assumed in the hydrodynamic calculations from 
which mass-loss formulas were derived. On the other hand, 
the uncertainties resulting from such inconsistencies are no 
more important than those resulting from a fundamental 



22 



Marigo et al.: Optical to far-IR isochrones 



lack of understanding of certain aspects of stellar evolution tensive testing of our TP-AGB tracks and dust models against 
and mass loss. observational data. 



Work is in progress to verify or improve on the points above; 
in particular, we are performing an extensive study of the 
Magellanic Cloud data (from / to [24], in preparation), and 
checking the results at low metallicities using data for dwarf 
galaxies and globular clusters (Gullieuszik et al. 2007ab, and 
work in preparation). Model atmospheres of cool giants are be- 
ing extended by Aringer et al. (in prep.). 

Despite the caveats, we believe that our new isochrones 
constitute a significant progress with respect to the past. We 
recall that the simple discussion of the above mentioned points 
would have been meaningless in the context of the more crude 
TP-AGB isochrones listed in Table 1 . 

In fact, we believe this is the first extended set of isochrones 
that includes the TP-AGB evolutionary phase with the level 
of completeness and detail it deserves. Crucial effects such 
as the third dredge up (which makes C-type stars), hot bot- 
tom burning (which largely prevents the most massive AGB 
stars from becoming C-type), and variable molecular opacities 
(which causes the cool red tail of C stars in near- and mid-IR 
CMDs) have been considered. Mass-loss rates are described 
with a formaUsm that discriminates between C- and M-type 
stars, and leads to the attainment of the super- wind regime dur- 
ing the latest stages of the TP-AGB evolution. The transition 
from the first overtone to the fundamental mode is predicted, 
which helps triggering higher mass-loss rates. Parameters asso- 
ciated with the third dredge-up are calibrated so as to reproduce 
the C-star luminosity functions observed in the Magellanic 
Clouds, and the lifetimes of C- and M-type stars derived from 
star counts in Magellanic Cloud clusters. Many photometric 
systems are considered, properly taking into account the dif- 
ferent spectra of 0-rich and C-rich stars, including the stellar 
obscuration by circumstellar dust at optical to near-IR wave- 
lengths, and its dust emission in the mid- and far-IR. The ex- 
pansion velocity of the wind and the dust-to-gas ratio are made 
to vary consistently as a function of stellar parameters (M, Teff , 
L, C/O, Z), in place of simpler parametrizations often used in 
the literature. Dust reprocessing of radiation makes the most 
evolved AGB stars to disappear in optical pass-bands, and ap- 
pear at a very extended red tail in near- and mid-IR CMDs. All 
these features make the new isochrones potentially very useful 
in the study of systems containing AGB stars, both resolved 
and unresolved. 

The new models have been included in the TRILEGAL 
code (Girardi et al. 2005; Girardi & Marigo 2007blil for the 
simulation of resolved stellar populations in the Milky Way and 
in external galaxies, and in the GRASIL code (GRAphite and 
SILicate; Silva et al. 1998, Panuzzo et al. 2005l3for comput- 
ing the integrated spectrum of galaxies from the X-ray to the ra- 
dio domain taking into account the effects of dust. Forthcoming 
papers will present applications of both codes, including the ex- 



'^ http : //trilegal . ster . kuleuven . be/J 
or http://st ev . oapd . i naf.it/trilegal 

'* |http : //a dlibitum . oat . ts . astro . it/ silva/def ault . html 
or |http : //web . oapd . inaf . it/granato/grasil/grasil . html | 
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Appendix A: Revised evolution of massive 
TP-AGB models 

The TP-AGB evolution of more massive stars, with initial 
masses 3.0 < M/Mq < 5.0, has been modified as follows. The 
starting point that has stimulated the revision is the behaviour 
of the integrated colours (e.g. V - K, J - K, H - K) as a func- 
tion of age as displayed by the Magellanic Clouds' clusters of 
age t k, 10** yr As shown in Fig. lA.ll the comparison with SSP 
colours, derived from the isochrones including the original TP- 
AGB tracks from Paper I, would indicate a sizable excess in 
the predicted V - K colours at ages around 1 - 4 x 10^ yr. The 
sudden jump to V - A" ~ 2.5 at f ~ 10** yr is due to the oc- 
currence of the AGB phase in models with Mi ~ 5 Mq, which 
corresponds to the maximum initial mass (for our set of mod- 
els) required to develop an electron-degenerate C-O core. In the 
1 -4 X 10^ yr age range, while the integrated V-band luminosity 
is essentially produced by less evolved stars still populating the 
region close to the main-sequence turn-off (see e.g. Chariot & 
Bruzual 1991), a significant contribution in the /T-band comes 
from the most massive AGB stars, with initial masses in the 
range 3.0Mq < M; < 5.0Mq. 




log(t/yr) 

Fig.A.l. Integrated V - K colours as a function of SSP age, 
derived from the original Paper I TP-AGB tracks (dashed line), 
and after the revision applied to the TP-AGB models with ini- 
tial masses > 3.0 Mq (solid line). The adopted metallicity is 
Z = 0.008. Observed colours of LMC clusters are taken from 
the compilations by Persson et al. (1983; filled circles), Kyeong 
et al. (2003; empty squares), and Pessev et al. (2006, starred 
symbols). 



It follows that the V - K colour jump at f ~ 10 yr should 
be ascribed to an overestimation of the TP-AGB phase of more 
massive AGB models, that translates into either a) too high lu- 
minosities and/or b) too long lifetimes. Of course these aspects 
are inter-related as they are both highly sensitive to model de- 
tails dealing with two crucial processes that affect the evolution 
of these stars, namely HBB and mass loss by stellar winds. 

Given the high degree of uncertainty that affects the the- 
oretical predictions for M along the AGB phase, the theoreti- 
cal discrepancy is likely attributed to an underestimation of the 
mass-loss efficiency in massive AGB stars. A simple remedy 
could be obtained artificially just introducing a multiplicative 
factor in front of our mass-loss law (similarly to the rj param- 
eter of the Reimers law) for Mi ^ 4.0 M©, but in this case we 
would introduce a discontinuity in our adopted formalism for 
M, as the correction factor should then vanish at lower Mi. 

After a general re-consideration of the many prescriptions 
of our synthetic TP-AGB model, we have found that a more 
physically-sound way to obtain higher mass-loss rates may 
come from an improved treatment of the pulsation properties 
of the most massive AGB stars, compared to the original de- 
scription given in Paper I. 

Specifically, following the results of non-adiabatic pulsa- 
tion models of Fox & Wood (1982)for stellar masses 3-6Mq, 
it seems likely that in Paper I the switching of the pulsation 
mode from the first overtone to the fundamental one has been 
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Fig.A.2. Transition luminosity defining the switching from 
fundamental mode to first-overtone pulsation mode, as a func- 
tion of the current stellar mass, assuming Z = 0.008 and 
log Teff = 3.52. The original relation of Paper 1 is shown as a 
solid line (Eq. lA.lb . while the modified relation at higher stellar 
masses (Eq. lA.2l l is drawn with a dashed line. The starred sym- 
bol corresponds to a pulsation model by Fox & Wood (1982) 
with Mi = 3.0 Mo, Z = 0.019, log(reff) = 3.49), while the ver- 
tical bar roughly defines the range of critical luminosity for a 
few models calculated Fox & Wood (1982) with Mi = 6 Mq. 

assumed to occur at too high luminosities in the most massive 
TP-AGB models. 

We recall that the critical luminosity, Li_o, is defined as the 
luminosity at which the growth rates of the first overtone mode 
and the fundamental mode are equal. In Paper I the critical lu- 
minosity has been derived from the pulsation models calculated 
by OstHe & Cox (1986) in the mass range 0.8 - 2.0 Mg and the 
fitting relation reads 



log(Li_()/L0) 



MG07 



-14.516 H- 2.277 log(M/Mo) (A.l) 

H-5.046 log Teff - 0.084 log(Z/0.02) . 

In Figure [A2] we compare Li_() as predicted by the original 
Paper I prescription (solid line), extended over the entire mass 
range of interest, with the results from Fox & Wood (1982). We 
see that while at M = 3 M© Eq. ( lA.ll i is fully consistent with the 
results from Fox & Wood (1982), at M = 6 Mq, (Li_o/i'0)MGO7 
exceeds significantly the range of critical luminosities derived 
from the Fox & Wood (1982) data. Based on these results we 
apply a rough correction A log Lcm to Eq. dA.ll ) in order to de- 
crease the transition luminosity for M ^ 3.0 M©, so that we 
write 



The revised relation for Li-o/L© is also plotted in Fig. IA.2I 
(dashed line). We see that it coincides with the old relation for 
M < 2.5 Mq, while it flattens to lower luminosities for larger 
masses. It should be observed that, given the paucity of the the- 
oretical pulsation models presently available, this revision is 
meant to be just a crude attempt to improve our treatment of 
Li_o- Indeed, new pulsation models over wider ranges of the 
key parameters (i.e. mass, metallicity, C/O ratio, etc.) are ur- 
gently required to allow for a real progress in modelling the 
pulsation properties of AGE stars. 

The new formula for the critical luminosity Li_() is used 
to re-calculate the TP-AGB evolution of models with initial 
masses 3.0 < MJMq < 5.0 for all choices of the metallicity. 
Massive models start to pulsate in the fundamental mode at 
lower luminosities than assumed in Paper I, which favours the 
attainment of larger mass-loss rates with consequent earlier ter- 
mination of the AGE phase. This effect is more pronounced at 
larger masses (M > 4 Mq) shortening the TP-AGE lifetimes by 
a factor of 2 - 5. The impact on the integrated Y -K colours is 
good as we see in Fig. lA.ll The jump at f ~ 10** yr is reduced 
and a more gradual increase of the colour proceeds at larger 
ages, allowing for the seeked consistency with the observa- 
tions. Similar results are obtained for other integrated colours, 
like 7 - /:, and // - K. 

It should also be emphasized that the TP-AGE fuel at ages 
t ^ 10'^ yr cannot be arbitrarily cut down since otherwise we 
would run into other problems, related to the observed prop- 
erties of the most luminous AGE stars and their progeny. For 
instance, we have verified that even invoking an extreme short- 
ening of the TP-AGE phase for these stars the predicted in- 
tegrated Y - K and J - K colours are still consistent with 
the observed ones, but in this case we cannot fulfill the nu- 
cleosynthesis constraint set by type I planetary nebulae. These 
latter are characterized by high He and N abundances, which 
are currently interpreted as the result of hot-bottom burning in 
more massive AGE stars (e.g. Marigo et al. 2003b; Pottasch & 
Eernard-Salas 2006). A too short activation of the CNO cycle 
at the base of the convective envelope would not allow to reach 
the He and N surface enrichment measured in type I planetary 
nebulae. In other words, the nuclear fuel burnt during the TP- 
AGE phase cannot be constrained only on the base of the in- 
tegrated luminosities, but it should also be checked against in- 
dications provided by abundance analyses in the stellar ejecta. 
In fact, as extensively discussed by Marigo & Girardi (2001) 
a precise fraction of the nuclear fuel consumption (emitted 
light) in post-main sequence stars ends in the form of chemi- 
cal (mainly He,N,0) yields. One cannot produce large amounts 
of He and N in the most massive AGE stars (as observed in 
their PN descendants), without adding a marked contribution 
of these same stars to the integrated fight. 



log(Li_()/L0) = log (Li_()/Lo)mgo7 + AlogLc, 
with 



(A.2) 



A log Lc. 



-0.125 (M/Mo- 2.5) 



(A.3) 



